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ABSTRACT 

A robust analysis of galaxy structural parameters, based on the modeling of 
bulge and disk brightnesses in the BVRH bandpasses, is presented for 121 face-on 
and moderately inclined late-type spirals. Each surface brightness (SB) profile is 
decomposed into a sum of a generalized Sersic bulge and an exponential disk. The 
reliability and limitations of our bulge-to-disk (B/D) decompositions are tested 
with extensive simulations of galaxy brightness profiles (ID) and images (2D). 
We have used repeat observations to test the consistency of our decompositions. 
The average systematic model errors are ^20% and ^5% for the bulge and disk 
components, respectively. The final set of galaxy parameters is studied for vari- 
ations and correlations in the context of profile type differences and wavelength 
dependences. 

Galaxy types are divided into three classes according to their SB profile 
shapes; Freeman Type-I and Type-II, and a third "Transition" class for galax- 
ies whose profiles change from Type-II in the optical to Type-I in the infrared. 
Roughly 43%, 44%, and 13% of Type I, II, and Transition galaxies respectively 
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comprise our sample. Only Type-I galaxies, with their fully exponential disks, 
are adequately modeled by our 2-component decompositions and our main results 
focus on these profiles. We discuss possible interpretations of Freeman Type-II 
profiles. 

The Sersic bulge shape parameter for nearby Type-1 latc-type spirals shows 
a range between n =0.1-2 but, on average, the underlying surface density profile 
for the bulge and disk of these galaxies is adequately described by a double- 
exponential distribution. The distribution of disk scale lengths shows a decreasing 
trend with increasing wavelength, consistent with a higher concentration of old 
stars or dust (or both) in the central regions relative to the outer disk. We 
confirm a coupling between the bulge and disk with a scale length ratio {ve/h) = 
0.22 ±0.09, or (hbuige/hdisk) = 0.13 ±0.06 for late-type spirals, in agreement with 
recent N-body simulations of disk formation. This ratio increases from ~ 0.2 for 
late-type spirals to ~ 0.24 for earlier types. These observations are consistent 
with bulges of latc-type spiral galaxies being more deeply embedded in their 
host disk than earlier-type bulges, as discussed by Graham (2001). Bulges and 
disks can thus preserve a nearly constant r^/ h but show a great range of surface 
brightness for any given effective radius. The similar scaling relation for early 
and late-type spirals suggests comparable formation and/or evolution scenarios 
for disk galaxies of all Hubble types. In the spirit of Courteau, de Jong, & Broeils 
(1996) but using our new, more extensive data base, we interpret this result as 
further evidence for regulated bulge formation by redistribution of disk material 
to the galaxy center, in agreement with models of secular evolution of the disk. 

Subject headings: galaxies: spiral — galaxies: photometry — galaxies: structure 
— galaxies: formation — galaxies: simulations 



1. Introduction 

Stellar density distributions provide important constraints for bulge and disk formation 
models. Historically, astronomers have embraced the r^/^ brightness "law" (de Vaucouleurs 
1948) and exponential brightness profile^ (de Vaucouleurs 1959a; Freeman 1970) to model the 



^The exponential nature of galaxy disk profiles emerges naturally in analytical models of disk formation 
(e.g. Lin & Pringle 1987; Dalcanton, Spergel, & Summers 1997; Ferguson & Clarke 2001). 
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light distribution of the galaxy bulge and disk, respectively^. Departures from the standard 
de Vaucoulcurs profile in the inner light distribution of early- and late-type spirals have 
however been demonstrated in a number of early studies (de Vaucouleurs 1959; van Houten 
1961; Burstein 1979), including the Milky Way (Kent, Dame, & Fazio 1991). Andredakis & 
Sanders (1994), de Jong (1996a), Courteau, de Jong, & Broeils (1996), and Carollo (1999) 
later used small samples of high-quality surface brightness (SB) profiles to establish the 
exponential profile as a better match to late-type disk bulges; thus SB profiles of most 
late-type spirals are best modeled by a double-exponential fit to the bulge and disk. 

A broader analysis suggests a range of bulge shapes from early- to late-type spirals 
(Andredakis et al. 1995; de Jong 1996a; Courteau et al. 1996; Graham 2001). Most of these 
analyses rely on the modeling of a generalized surface density function such as that proposed 
by Sersic (1968); 



l/n 



or, in magnitudes. 



/(r) = /oexp { - 



/x(r)=/xo + 2.51og(e)<{ (-) }. (2) 



where Hq (Jq) is the central surface brightness (intensity), Tq is a scaling radius, and the 
exponent l/n is a shape parameter that describes the amount of curvature in the profile. 
For n = 1 or 4 one recovers a pure exponential or the de Vaucouleurs r^/^ profile respectively. 

Collectively, the works above suggest that the bulge shape parameter n correlates with 
absolute luminosity and half-light radius, such that bigger, brighter systems have larger 
values of n. This result was extended to brightest cluster galaxies by Graham et al. (1996). 
Courteau et al. (1996) also demonstrated a tight correlation between the bulge and disk 
exponential scale lengths, for all spiral types, with hb/hd — 0.1 ± 0.05 (where h = tq and 
n = 1 in Eq. 1). The exponential nature of late-type galaxy bulges and the correlation 
between bulge and disk scale lengths was interpreted by Courteau et al. (1996) as evidence 
for regulated bulge formation by redistribution of disk material to the galaxy center by a bar- 
like perturbation. We will return to this important constraint for secular evolution models 
in 5 5.3. 



^It is important to remind ourselves from the onset that bulge-to-disk decompositions, and inward ex- 
trapolations of the disk into the central bulge and/or bar, may have no physical (or dynamical) basis. They 
provide a convenient description of the light distribution of a galaxy's components that arc otherwise dy- 
namically coupled. The effective integrals of motions arc likely similar for all the co-spatial components, 
though kinematically distinct bulges (counter-rotating nuclei) are known to exist. 
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This study focuses on the development of a reliable set of observables and constraints 
for structure formation models. An important goal is to measure the range of the Scrsic 
n parameter for virialized disk systems. The analyses described above arc reproduced and 
expanded upon with the largest multi-band survey of its kind to date and a clearer under- 
standing of model limitations than previously attained. We aim to characterize and quantify 
the intrinsic structural properties of the bulge and disk and the extent of their variation 
with wavelength. These characterizations are made through reliable modeling of bulge and 
disk parameters from SB profile decompositions. Multi-wavelength information also pro- 
vides insight about structural variations within and among galaxies due to dust and stellar 
population effects. While some of these issues have been addressed before, there remains 
a number of significant measurement uncertainties and technical limitations which we now 
investigate thoroughly. 

This paper is organized as follows: a brief description of the database is given in § 2 
and in § 3 we discuss our B/D decomposition algorithms (ID and 2D) and the simulations 
to test the reliability of our technique. For the readers interested mostly in final profile 
decompositions and results, a summary of the simulation results and guidelines is given in 
§ 3.4. Actual B/D decompositions of galaxy SB profiles are presented in § 4, followed by a 
discussion and interpretation of the results in terms of secular evolution models in § 5. A 
discussion on the nature of Freeman Type-II profiles is also presented in § 5. We conclude 
with future directions in § 6. Two appendices present (A) a discussion of the functional form 
for the Sersic coefficient and (B) decomposition results for our Type-I profiles. 



2. The Data 

Our structural analysis of galaxy luminosity profiles is based on the catalog of multi-band 
images of late-type spiral galaxies by Courteau, Holtzman, & MacArthur (2002; hereafter 
Paper II). It consists of over 1000 deep B, V, R, and H images of 322 nearby bright late- 
type spiral galaxies. The data were collected between 1992 and 1996 at Lowell Observatory 
and Kitt Peak National Observatory (KPNO). A full description of the sample selection, 
observations, and reductions is presented in Paper II. A summary is given below. 

The galaxy sample was selected from the Uppsala General Catalogue (UGC, Nilson 
1973) with the following criteria: 

• Predominently late Hubble types 

• Zwicky magnitude < 15.5 

• Blue Galactic extinction = 4 x E{B - F) < 0^5 (Burstein & Heiles 1984) 
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• Inclination bins covering face-on {i < 6°), intermediate (50° < i < 60°), and edge-on 
{i > 78°) projections 

• Blue major axis < 3' . 

This catalog is not complete in any sense of the (much abused) term. The diameter limit 
was constrained primarily by the field of view of the infrared cameras in use at KPNO (IRIM 
and COB) and Lowell Observatory (OSIRIS) in 1992-1996 and the requirement for blank 
areas in the field of view for sky subtraction. Additionally, peculiar and interacting galaxies 
(e.g. no visible tidal tails) were excluded to ensure that the sample consisted only of isolated 
disk dominated galaxies. Barred galaxies, as classified in the UGC, were not excluded per 
se but only a handful were observed. For the present analysis, we use a sub-sample of 121 
galaxies with face-on and intermediate inclinations only, for a total of 523 images^. The 
distribution of Hubble types in our reduced sample is: 2 Sab, 26 Sb, 19 Sbc, 38 Sc, 25 Scd, 
11 Sd. 

All distances are corrected to the reference frame of the Local Standard of Rest (Courteau 
& van den Bergh 1999), and we use Hq = 70 km sec"^ Mpc"^. The survey effective depth is 
(cz) ~ 5500 km s"^ or 80 Mpc. 

2.1. Observations and Basic Reductions 

All optical BVR images were obtained from 1992 to 1994 at Lowell Observatory with 
a TI 800x800 chip (scale = 0'/5/pix) on the Perkins 72" telescope. The infrared H-band 
images were acquired from 1993 to 1995 at KPNO with the 2-meter and 4-meter telescopes 
equipped with either a HgCdTe (IRlM) or an InSb (COB) 256x256 array (l'.'09/pix and 
O'.'5/pix respectively), and from 1995 to 1996 with the OSIRIS imager (l'.'49/pix) mounted on 
the Lowell 72" telescope. The exposure times were: 300s in R, 400s in V, 1500s in B, and on- 
target integration of 1200s in H. Landolt (1992) standards covering a wide range of airmasses 
and colors were observed each night at Lowell Observatory, giving a photometric accuracy 
of ~2% for the optical passbands. UKIRT standards (Guarnieri, Dixon, & Longmore 1991) 
observed each night yielded H-band photometric calibrations good to ~3%. Stars and defects 
were edited from the images prior to further analysis. 

The typical seeing full-width at half-maximum (FWHM) at Lowell and KPNO was 
2"0 with typical standard deviations of ~20% (optical) and ~35% (IR) per image. These 



^Radial brightness profiles cannot be measured for fully edge-on galaxies. Our images of edge-on systems 
will be used in a forthcoming analysis of stellar and dust scale heights and truncation radii in spiral disks. 
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measurements were computed as the mean of the FWHMs of all non-saturated stars measured 
automatically on each image frame; typically 5 to 40 measurements per frame were used. 

We measure mean sky levels (for a 5-6 day old moon) oi B — 21.9 ± 0.8, V — 21.2 ± 
0.5, R — 20.6 ± 0.5, and H — 14.1 ± 1.2 mag arcsec"^ . Typical systematic errors in the sky 
measurement, computed from 4 or 5 sky boxes suitably located between the galaxy and the 
edge of the frame, are 0.5 — 1.0% in the optical and 0.005 — 0.01% in the JR. 

Azimuthally-averaged SB profiles were extracted for all the galaxies using ellipse fitting 
with a fixed center. To ensure a homogeneous computation of structural parameters and 
color gradients, we use the isophotal maps from the R-band to determine the SB profiles 
in BVH. Even though dust effects can still play a role at 7000 A, the R-band was adopted 
for our isophotal templates as it has the most stable sky and deepest profiles. Wc allowed 
a variable position angle and ellipticity at each isophote, but a comparison with SB profiles 
extracted using concentric isophotal fits demonstrated that our results do not depend on the 
fitting technique. Further information about profile extraction and CCD surface photometry 
can be found in Courteau (1996a) and Paper II. We trace SB profiles to ~26 mag arcsec"^ 
in optical bands and ~22 mag arcsec"^ at H-band. These levels correspond to a surface 
brightness error of ~0.12 mag arcsec"^ . 



2.2. Surface Brightness Corrections 

The observed surface brightness of a galaxy can change when viewed at different inclina- 
tion angles, depending on the distribution of a galaxy's interstellar medium and its opacity. 
Surface brightnesses are also affected by Galactic foreground extinction and redshift dim- 
ming. We account for the latter effects but defer any treatment of internal extinction, which 
vary greatly from author to author, to Paper II. Our conclusions do not depend on the exact 
values of the central and effective surface brightnesses of galaxies. 

We correct for Galactic foreground extinction using the reddening values. Ax, of Schlegel 
et al. (1998) and assuming an Ry — 3.1 extinction curve (e.g. CardeUi et al. 1989), 

f^c,Gal = /^ofes ~" ^A- (3) 



We correct surface brightnesses for the (1 -|- cosmological redshift dimming [per unit 
frequency interval) as 

l^i, = l^obs-7.5\og{l + z). (4) 
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The final correction to the observed surface brightnesses is thus, 

/x^ = /xl-/lA-7.51og(l + ^). (5) 

Examples of the types and extent of the SB profiles for typical late- type spiral galaxies in 
our sample are shown in Fig. 1. For Type-I disks (Freeman 1970), the inner profile always lies 
above the surface brightness of the inward extrapolation of the outer disk, whereas Type-II 
systems have a portion of their brightness profiles lying below the inward disk extrapolation. 
We define a Transition case for luminosity profiles that change from Type-II at optical 
wavelengths to Type-I in the infrared. Many galaxies classified as Type-II show a weakening 
of the inner profile dip at longer wavelengths and, in this sense there is no clear distinction 
between the Type-11 and Transition galaxies. Likely interpretations for Type-II profiles are 
discussed in § refsubsec:typell. 

3. Simulations of Bulge-to-Disk Decompositions 

In order to measure galaxy structural parameters, we have developed two independent 
algorithms to decompose the galaxy ID and 2D light distributions into bulge and disk com- 
ponents. These programs allow for a generalized Sersic bulge, an exponential disk, and a 
central bar for 2D images. There are several issues involved with accurate decompositions, 
particularly with the measurement of bulge parameters, including; the sensitivity of final 
results to starting guesses, effects of statistical and systematic errors in sky brightness and 
seeing estimates, choice of fit basehne, etc. We explore these in great detail below using both 
ID and 2D analyses to determine the robustness of our codes and the reliability of our final 
solutions. Because projected surface brightness profiles contain fewer data points than full 
2D images, we can create ID simulations faster than 2D models. Thus our most extensive 
tests rely on ID simulations, which are shown to be fully consistent with 2D simulations 
when considering axisymmetric features. 

3.1. ID and 2D Algorithms 

Our brightness profile (ID) bulge-to-disk (B/D) decomposition algorithm was initially 
developed by Broeils & Courteau (1997) and subsequently improved by LM. This program 
reduces ID projected galaxy luminosity profiles into bulge and disk components simultane- 
ously using a non-linear Levenberg-Marquardt least-squares (NLLS; see § 15.5 in Press et 
al. (1992)) fit to the logarithmic intensities (i.e. magnitude units). Random SB errors are 
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accounted for in the (data— model) minimization, whereas systematic errors such as uncer- 
tainties in the sky background and determination of the image mean PSF are accounted for 
separately in a series of experiments designed to calibrate their effects. Seeing effects in our 
model galaxies are accounted for by convolving the theoretical bulge-disk surface brightness 
profiles and images with a radially symmetric Gaussian Point Spread Function (PSF), of the 
form 



where Itotaii^) is the intrinsic surface brightness profile, a is the dispersion of the Gaussian 
PSF and Jq is the zero-order modified Bessel function of the first kind (see also Trujillo et 
al. (2001) for a study of the Moffat PSF.) 

The 2D decomposition program is based on the same NLLS technique as above but 
uses the full 2D image in intensity units instead of a logarithmic radial surface brightness 
profile. While computationally more intensive than its ID analogue, the 2D decomposition 
should yield more physically meaningful results since the azimuthal information is lost in ID 
profiles. Byun & Freeman (1995), de Jong (1996a), and Simard et al. (2002) have discussed 
the merits of the 2D approach, such as greater ability to recover true parameters (based on 
simulations), and the potential to model non-axisymmetric features such as bars, rings, and 
spiral arms. The need for the implementation and testing of a robust 2D B /D decomposition 
package is thus obvious, but we find that ID decompositions compare favorably for reliability 
and predictive power provided high S/N ID radial profiles are used. Note that neither ID nor 
2D decompositions are impervious to dust extinction effects. Extinction effects are lessened 
at H-band, but can still be significant in disk bulges and spiral arms. A proper recovery of 
the true stellar density profile would require a full 3D radiation transfer treatment, and such 
an analysis is beyond the scope of this work. 



A fundamental aspect of profile decompositions is the choice of fitting functions. The 
disk fight is modeled with the usual exponential function. 





3.2. Methodology 




(7) 



or, in magnitudes. 




(8) 
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where //q = — 2.51og/o and h are the disk central surface brightness (CSB) and scale length 
respectively, and r is the galactocentric radius measured along the major axis. In the 2D 
decompositions, the computation of the radius at each pixel requires two additional param- 
eters: the position angle (PA) of the disk major axis on the sky and the disk ellipticity, 
e = 1 — b/a, where a and b are the major and minor axes of the disk respectively). To test 
for the shape of the bulge luminosity profiles we adopt the generalized formulation of Sersic 
(Eqs. 1 & 2). 

It has become customary to express the disk parameters in terms of scale length and 
CSB {h and ^q), while the bulge parameters arc expressed in terms of effective parameters 
(re and /Xg)- We adopt this formalism, thus parameters with subscript e refer to the bulge. 
Eq. 1 can be re-written as: 



Ib{r) = /ecxp 



1/n 



(9) 



where the effective radius, r^, encloses half the total extrapolated luminosity^, /g is the 
intensity at this radius and 6„ is chosen to ensure that 



roo pre 

/ Ib{r) 27rr dr ^ 2 h{r)2T:rdr. 
Jo Jo 



(10) 



In magnitudes Eq. 9 translates to 

libir) ^ lJ>e + 2.5 log(e) 6„ 



l/n 



- 1 



(11) 



where /ig is the effective surface brightness. 

It is trivial to convert from Eq. 2 to Eq. 11 by noting that 

^jbe = A^o + 2.5 log(e) bn- 



(12) 
(13) 



Eq. 10 implies that 

T{2n) = 2-1 {2n.br,) (14) 

where r(a) is the gamma function and 7(0, x) is the incomplete gamma function. Unfor- 
tunately, Eq. 14 cannot be solved analytically for 6„. Various numerical approximations 



^For a pure exponential disk, = 1.678/i. 
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have been been given in the hterature (Caon et al. 1993; Graham Sz Prieto 1999; Ciotti & 
Bertin 1999; Khosroshahi, Wadadekar, & Kembhavi 2000; MoUenhoff & Heidt 2001). One 
often encounters the approximation 6„ ^ 2n — 0.32, vahd supposedly for all values of n (sic). 
Khosroshahi et al. (2000) contend that this approximation is accurate to one part in 10^, 
with a range of validity on n unspecified. However, because the gamma function diverges 
near the origin, most utilized approximations are inaccurate for values of the Sersic expo- 
nent n < 1. Differences between numerical solutions for 6„ (Eq. 14) and commonly adopted 
approximations can yield brightness differences greater than 0.1 mag arcsec"^ for n < 1. 
As we wish to test for bulges with Sersic n parameter smaller than 1, we have adopted the 
asymptotic expansion of Ciotti & Bertin (1999) to 0{n~^) for n > 0.36. For n < 0.36 this 
solution diverges and instead we use a polynomial expression (4*^ order) accurate to one part 
in 10^. We compare different numerical solutions for 6„ (Fig. 23) and present our adopted 
functional form in Appendix A. 

An ilhistration of profile shapes for different values of the Sersic n parameter is shown 
in Fig. 2. The top panel shows profiles with He = 21 mag arcsec^^ and rg = 3'.'5 for values 
of n in the range 0.2 < n < 4. The bottom panel shows the same profiles but for a constant 
CSB of fiQ = 18 mag arcsec"^ . For n < 1 the profiles are shallow at small radii (< r^) and 
fall off rapidly with increasing radius. Conversely, profiles with n > 1 are steep at small 
radii (<C rg), but level off as r increases. Given the large differences in the profile shapes 
above and below n — 1 (exponential case), one might expect different physical mechanisms 
(formation, transport, dynamics, interactions) to be at work for systems whose light profiles 
have very different n values. Additionally, for the small bulges of late-type galaxies, poor 
seeing could conceivably smear the image such that an intrinsically n > 1 bulge could be 
mistaken for an n < 1 structure. 

Of potential relevance to the study of galaxy structure is the relative fight fraction 
contributed by the bulge and disk. This is expressed in terms of a bulge-to-disk luminosity 
ratio, B/D, derived by integrating the bulge and disk luminosity profiles to infinity. For a 
face-on Sersic profile the total extrapolated luminosity is given by 




(15) 



and for a face-on exponential disk 




(16) 



giving a bulge-to-disk fight ratio of 




(17) 
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Eqs. 15 & 16 should be multiplied by the factor {h/a) when considering projections on 
the plane of the sky. One may use Eq. 17 in a general sense, independent of projection, 
under the assumption that the bulge and disk density distributions have similar axes ratio 
(nearly true for late-type galaxies). A weakness oi B / D ratios for systematic comparisons 
of galaxy light profiles is its model dependence and the potential covariances between some 
of the model parameters. Consider Fig. 2 (top) for the relative light fractions contributed 
by profiles of different n values, normalized to n = 1. The integrated bulge light increases 
steadily as a function of n, for given values of and //g- Thus, the adopted n value in a bulge- 
to-disk decomposition has a strong infiuence on the computed B/ D ratio. Additionally, since 
larger n profiles contribute light out to large r, the combination of a high n and a low //g 
(bright r^/^ bulge) could take away light from the outer disk and artificially boost the B /D 
ratio. A discussion on non-parametric statistics, such as concentration indices (Kent 1985; 
Courteau 1996a; Graham 2001) which alleviate model dependences, is presented in Paper II. 

We model the total galaxy luminosity profile as a sum of bulge -|- disk components: 

Itot{r)^h{r) + h{r). (18) 

Profile smearing by atmospheric turbulence is accounted for in B/D decompositions by con- 
volving Eq. 18 with a Gaussian PSF of the form of Eq. 6. 

Similar B/D analyses have also considered additional terms for a Gaussian bar (de Jong 
1996a), a lens or ring (Pricto ct al. 2001), spiral arms, and stellar disks with inner and/or 
outer truncations (Kormendy 1997; Baggctt, Baggctt, & Anderson 1998). We restrict our 
choice of fitting functions to a Scrsic bulge and a non-truncated exponential disk for a 
number of reasons. We find no prominent bars in our sample and most our disk profiles are 
fairly linear (in magnitude space). Azimuthal averaging for ID profiles smoothes out spiral 
arm features (to a different extent depending on whether the position angle was fixed or 
allowed to vary in the profile extraction. Removal of spiral arm signatures from the light 
profiles or images would require more time and effort than is warranted by our analysis at 
this stage.) We do not consider a sharp inner disk truncation for a number of reasons: (i) 
unsharp masking techniques reveal spiral structure from the inner disk into the galaxy center 
(Courteau 1992, 1996b; Elmegreen, Elmegreen, & Eberwein 2001); (ii) using HST images 
of inner disks, CaroUo (1999) also finds evidence for inner spiral structure and nuclear star 
clusters in the centers of early- to intermediate-type spiral galaxies; (iii) all components of 
the Galaxy have their peak surface brightnesses in the center (e.g. Wyse 1999). Thus, at 
least some evidence suggests that spiral disks reach in all the way to the center of late-type 
systems. A lowering of the disk central surface density may occur as stars get heated up 
into a bulge by the action of a bar-like instability. An exponential profile with a core may 
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thus be a reasonable description of Type-II SB profiles. We do not consider this approach 
here, but point out that resolved B/D kinematics of nearby galaxies would provide a clear 
indication whether stellar populations have been strongly depleted and/or systematically 
scattered vertically into a bulge. Some of our galaxies show outer disk truncation (e.g. 
Fig. 13), but see § 4.3.2. 

The best-fit parameters of the (data— model) comparison are those which minimize the 
reduced chi-square merit function, described in intensity units as 



where N is the number of data points used. M is the mimber free parameters (i.e. N - M = 
V = Degrees of Freedom), and (jj is the statistical intensity error at each pixel (2D) or surface 
brightness level (ID). From here on the u subscript will be omitted and the variable refers 
to a psr degree of freedom (unless otherwise specified). 

The global of intensities is clearly dominated by the contribution from the disk, 
virtually irrespective of the fitted bulge. This effect would be accentuated in galaxies with 
prominent features, such as spiral arms, rings, or lenses, which are not accounted for in our 
pure exponential disk models. Cases are found where B/D decompositions with significantly 
different bulge exponent n values for a given profile have nearly the same global value (see 
Figs. 8 & 9). Thus, in order to refine our parameter search for the best-fit bulge and disk 
model, we compute a separate, inner, statistic out to twice the radius where the bulge 
and disk contribute equally to the total luminosity of the galaxy {rb=d = 2r(/(, = Id)). We 
label this statistic as Xin (see Graham 2001 for a similar formulation). For cases where the 
bulges are so small that they never truly dominate the light profile (i.e. rb=d is undefined), 
we compute the Xin ^^e radius at which v = 1. 



This section describes extensive testing of our bulge-to-disk decomposition programs. 
Artificial SB profiles and images were created with a wide range of bulge profile shapes 
and exponential disks including realistic noise and seeing effects. Real galaxies are clearly 
more complicated than the sum of two idealized mathematical functions, but these tests 
provide a reasonable base for a global understanding of the reliability and limitations of 
B/D decomposition algorithms. The mock catalog of SB profiles and images will be used to 
address the following questions: 




(19) 



3.3. Reliability of the Decomposition Results 
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• How reliable and meaningful are the bulge-to-disk decompositions and fitted parame- 
ters? 

• How crucial are initial estimates? Are model fits always converging to the lowest 
minimum? 

• How do seeing effects and sky subtraction errors affect the decompositions, and can 
they be properly accounted for? 

• Are the small bulges in late-type disk galaxies sufficiently resolved to permit a reliable 
solution of the Sersic n parameter as a free parameter? 

The literature abounds with investigations of profile fitting algorithms based on artifi- 
cial data, such as Schombert & Bothun (1987; hereafter SB87) who performed double-blind 
experiments where one of the authors created mock luminosity profiles and the other indepen- 
dently fitted the data. The SB profiles combined a de Vaucouleurs bulge and an exponential 
disk. Photon noise, at a level matching typical blue CCD performances, and a systematic 
0.5-3.0% error of the sky background were added to the profiles. SB87 found that the simul- 
taneous fitting of disk and bulge using standard NLLS techniques could reproduce the input 
parameters to within 10-20% in cases where galaxy profiles can be decomposed perfectly 
as the sum of a bulge and disk (which fails for Typc-11 profiles.) SB87 claim that a sky 
estimate uncertainty of up to 3% does not affect their derived parameters, but we find that 
sky errors as small as 1% can have a significant effect on the shape of the outer disk profile 
and the derived bulge and disk parameters (see § 3.3.5 below). SB87 did not consider other 
fitting functions but recognized that bulges may not be adequately described by the de Vau- 
couleurs r^/^ function. Andredakis & Sanders (1994) later examined the inadequacy of the 
T^l^ functional form for the bulge (ID) profile, and first established the double-exponential 
nature of late-type spirals. 

2D B/D decomposition techniques, which exploit the full galaxy image, were also de- 
veloped and tested in similar fashion in the mid-nineties (Byun & Freeman 1995; de Jong 
1996a). De Jong performed extensive tests with mock galaxies modeled as pure exponen- 
tial bulges and disks, exploring the effects of errors in the measured observables including 
the seeing FWHM, sky background level, minor over major axis ratio, hja. and position 
angle, PA. These observables were used as fixed input parameters to the fitting routine, 
and de Jong calibrated the effect of measurement error on the determined parameters by 
decomposing the artificial galaxies using erroneous values for each observable. He concluded 
that: errors in /xq are predominantly caused by sky subtraction errors and can be as large 
as 0.1 mag arcsec"^ ; errors in h can reach 10% and are dominated by sky background and 
ellipticity measurement errors; bulge parameter errors, of order 20%, are controlled by the 
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B/D size and brightness ratios. Bright bulges are most affected by seeing errors, and fainter 
bulges can also be affected by sky background errors. 

Our own investigation reaches similar conclusions and further extends de Jong's sim- 
ulations. We test for the robustness of the fitting procedure and accuracy of the derived 
parameters with various values of the fit initial estimates, seeing FWHM, sky value, and 
their errors, and - unlike de Jong (1996a) - we model the bulge with a generalized Scrsic 
profile. These simulations were initiated by Broeils & Courteau (1997) but are extended 
here in much greater detail, especially with respect to the determination of the bulge shape 
parameter n. 



3.3.1. Simulated Profiles and Images 

Our tests use a large set of artificial SB profiles and images which span a wide range 
of the bulge, disk, and seeing parameters. The mathematical forms of the bulge and disk 
components are those discussed in § 3.1. Noise was added to the model profiles and images 
from a Gaussian distribution with deviation representative of the standard brightness errors 
of our luminosity profiles at a given surface brightness level (see Courteau 1996a, Fig. 9; 
Paper II) . 

One hundred SB profiles and fourty images with realistic noise were created for each 
bulge, disk, and seeing combination. Most of the simulated profiles and images had the same 
disk parameters, 

jiQ — 20 mag arcsec"^ , 
h = 12", 

which are representative of a typical galaxy in our sample^. For exponential {n — 1) bulges, 
the structural parameters were selected from 

He — 16, 17,..., 22 mag arcsec"^ , 
re = 0.1, 0.2,..., S'.'O 

corresponding to B/D ratios ranging from to 5 and B /T from to 0.8 (see Eq. 17). 
All profiles and images were convolved with a seeing disk of 



^Different values of h were also tested but the results do not change significantly. 
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FWHM = 1.0, 1.5,..., 3'/0. 

Our models span the full range of parameters typically found in late-type bulges (e.g. de 
Jong (1996b); Courtcau (1996a); Broeils & Courteau (1997)) and the seeing values match 
the expected range at Lowell Observatory and KPNO (see § 2). 

We also explored the following range of Sersic n values 

n = 0.2, 0.4,..., 4.0 

for the set of combinations with = 0.8, 1.5, and 2'.'5, fie = 18, 20, and 22 mag arcsec"^ 
and for seeing FWHMs of 1.5, 2.0 and 2'.'5. Here the corresponding B/D ratios range from 
about to 1, and B/T from to 0.5. We also simulated the full range of Tg, FWHM, and 
jie = 18, 20, 22, and 24 mag arcsec"^ for n — 0.2 and n = 4.0 for a test regarding the initial 
estimates (see § 4.2). 

A total of about 223,000 artificial SB profiles and 16,200 images were created and 
modeled. The mock profiles were all sampled at O'.'5/pixel to match the data. 

The ID and 2D decompositions of the artificial profiles and images were deemed satis- 
factory if they met the following (rather liberal) criteria: 

• solution found within 100 iterations 

• < re(/it) < 500" 

• < He{fit) < 30 mag arcsec"^ 

• 0.05 < n^fit) < 10 

• < h^fu) < 100" 

• < ijLo{fit) < 30 mag arcsec"^ 

For each set of parameters, we require that at least two-thirds of the 100(40) simulated 
profile (image) decompositions pass these criteria to be included in the analysis. The 2D tests 
were not developed as fully due to prohibitive computing times. Thus our tests rely more 
heavily on the ID technique, but we have confirmed that the results from both techniques 
corroborate each other. 
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3.3.2. Disk Initial Estimates 

NLLS algorithms require initial estimates as input parameters, and we must verify 
whether our final solutions are sensitive to our initial guesses. Different initial estimates 
may yield different solutions with comparable values especially if the topology of the 
distribution is non-trivial or shallow (e.g. Schombert & Bothun 1987; de Jong 1996a). Our 
simulations confirm the robustness of our algorithms to a wide range of initial disk parameter 
estimates. The results are shghtly more sensitive for large values of n, but in general the 
disk parameters were perfectly recovered independent of the initial guesses. One must still 
caution that if the bulge is fit with the wrong n, the fitted disk parameters will differ from 
their intrinsic values, even if the initial estimates were good. Fig. 3 shows the relative fit 
error for /i, A/i, where 

^ /i-fit(niean) - /imodei ^2o) 

^model 

versus the fitted n value (held as a fixed parameter in the decomposition) for bulges with n 
— 0.2, 1.0, and 4.0. (i.e. A/i = indicates a perfect recovery of the model parameter h). 

In these decompositions the correct initial values were given for the bulge and disk 
parameters and the seeing FWHM was assumed to be known. Different seeing values are 
represented by three point types: circles, triangles, and squares for seeing values of 1.5, 2.0, 
and 2'/5 respectively. We see that the fitted disk parameters can have errors as large as 10% 
when the bulge is modeled with an incorrect shape parameter. The reason for this is obvious 
if one considers the different shapes for the different values of n shown above in Fig. 2, where 
the different shapes contribute differently to the outer profile. The bulge is even more 
sensitive to the fitted n (relative errors in excess of 300% for the worst cases, the figure is 
not shown). 



3.3.3. Bulge Initial Estimates 

We now investigate the importance of the bulge initial estimates. Given that the fit 
baseline for late-type bulges is much smaller than that of the disk and the apparent size of 
these bulges is comparable to the seeing disk, the bulge model is likely to be much more 
sensitive to the input parameters. The sensitivity to the initial estimates of bulge parameters 
re and /Xg was tested using offsets of ± 25% and ± 50% from the model, as well as the correct 
model values as initial estimates. This range of offsets matches typical estimate excersions 
in our data modeling (see § 4.2). The Sersic n was presumed known and the tests were 
performed for three different values of n =0.2, 1, and 4. 
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The tests for the n — 1 case demonstrated that bulge parameter initial estimates are 
not important in both the ID and 2D decompositions, as long as Vg > 0.3*FWHM. Below 
this limit initial estimates are more important; errors on r^. can exceed 50% and errors of up 
to A/Xg ± 0.3 mag arcsec"^ can occur. In the n = 0.2 case the parameter recovery is largely 
independent of the seeing FWHM, as expected for profiles that arc flat in the center. Given 
incorrect initial guesses, recovered parameters for profiles with > I'.'O can still be trusted, 
but for smaller bulges the algorithm is trapped in a local minimum and the output value is 
nearly the same as the input, e.g. a ±25(50)% input error yields a ±25(50)% output error. 
For n — 4.0 profiles, the parameter recovery is strongly dependent on the seeing FWHM 
such that decomposition results for profiles with rg < 0.7 * FWHM cannot be trusted. Here 
again input and output errors are approximately equal. Moreover, even with correct initial 
estimates, the model parameters are not perfectly recovered for profiles with < I'.'O. 

We have interpolated these results for different values of n and define a parameter space 
for which our solutions are not affected by the choice of initial estimates: 



The corresponding results for the 2D decomposition algorithm closely match those from 
the ID tests. 



The effect of an uncertainty in the seeing FWHM measurement on the model parameters 
(n, Te, yUe, h, and iiq) must be accounted for in our simulations. We fit each model using 
not only the nominal seeing value, but we also varied the seeing FWHM by typical seeing 
measurement errors (1 a; ~ 15 — 20% for optical, ~ 35% for infrared; see Paper 11). Figs. 4 & 
5 show the effect of an incorrect seeing estimate on the fitted for the ID and 2D algorithms 
respectively. Plotted are the relative fit error on rg, Avg, where 



re > (0.3)^/" * FWHM 



and 




-0.75 n + 1.15 for n< 1.0 
0.2 n + 0.2 for n> 1.0 



3.3.4. Seeing Effects 



Are = 



re,fit(mean) - re,modei 



(21) 



^e, model 



The ID and 2D simulations agree very well. Figs. 4 & 5 show that when the correct seeing 
is used as input, the bulge and disk parameters are recovered perfectly for the full range of 
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bulge parameters and for all values of the seeing FWHM tested. However, even moderate 
seeing uncertainties can severely affect bulge parameters depending on the size of the bulge 
relative to the seeing FWHM. If the seeing width is under(over)-estimated r^. is systematically 
over(under)-estimated, worsening for smaller and fainter bulges and larger seeing values. 
Similar trends are seen for fie- Our tests show that the fit errors can be significantly larger if 
the seeing FWHM is over-estimated than if it is under-estimated. As a rough rule of thumb, 
for Tg ~ FWHM and a seeing measurement uncertainty at the 35% (15%) level, the bulge 
can be trusted to within 10-25% (0-10%), and //g to within ± 0.1-0.4 (0-0.2) mag arcsec"^ , 
the lower end of the range applying to the brightest bulges and increasing towards the upper 
end for the fainter bulges. For rg ~ FWHM -|- 1 the errors improve to within 0-15% (0-10%) 
for Te, and ± 0.0-0.2 (0.0-0.05) mag arcsec"^ for He- 
There is no appreciable effect due to seeing on the disk parameters (less than 1%) except 
for the worst case of a FWHM of 3'.'0 and a 35% seeing over-estimate. In all other cases, the 
disk parameters are virtually unaffected by seeing, as the size of the disk is much larger than 
the seeing profile. However, it is of paramount importance to use accurate seeing estimates 
and realistic seeing errors in order to sample the true range of bulge parameters in B/D 
decompositions. 

3.3.5. Sky Uncertainty Effects 

We now test for the effects of an improper sky subtraction on the decompositions. The 
tremendous sensitivity of B/D decomposition and scale length determinations to sky errors 
has been highlighted before (Courteau 1992; de Jong 1996a). Here we aim to provide a firm 
quantitative assessment of such errors. We re-model the same simulated profiles as in the 
previous section but using sky values that are ±1% of the nominal sky level (typical error in 
the optical passbands), and using a typical optical sky brightness of 21 mag arcsec"^ ^. Since 
bulge brightnesses are typically greater than the sky level, at least at optical wavelengths, 
one might expect bulge parameters to be somewhat insensitive to sky subtraction errors. 
However, the outer disk is very sensitive to sky subtraction errors and a modified disk 
ultimately affects bulge structiuc due to their coupling. Quantitatively, our tests show that 
if the sky is over- or under-estimated by 1%, the error on the disk scale length, A/i, will be 
of order 5-15% and the disk CSB, A/xq, will be ± 0.1-0.25 mag arcsec^^ . These dispersions 
hold for the full range of bulge brightnesses except the two faintest bulges which are one and 



^Note that the percent sky error in the H-band is much smaher than in the optical, but the sky at H is 
much brighter, so the effect should be comparable. 
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two magnitudes fainter than the sky; these had errors in excess of 50%. The error ranges 
are controlled by the relative sizes of the bulge and disk such that the disk parameter errors 
increase slightly from smaller to larger bulges. This is simply because a larger bulge weakens 
the importance of the disk in the central parts, thus giving more weight to the sky-sensitive 
outer disk. 

The errors on the bulge parameters are negligible for bulges with (/i.,fcy — fJ-e) > 1 
mag arcsec"^ for the entire range of Tg and seeing FWHMs, but increase up to Ar^. > 15% 
and |A/ie| > 0.1 mag arcsec^^ (increasing as the bulge gets smaller and as seeing conditions 
degrade) for bulges with {fJ^sky — fJ'e) < 1 mag arcscc"^ . In other words, if the bulge effective 
surface brightness is less than one magnitude greater than the sky brightness, the bulge 
parameters will be strongly affected by sky subtraction errors. This effect is often neglected 
in studies of bulge/disk structure. 

The bulge and disk parameters are most affected for the case of an under-subtracted 
sky. This is largely due to a magnitude threshold of 26.5 mag arcsec"^ in our decomposition 
algorithm. The data are too noisy below this value (Paper II) and we exclude them from 
the fits. This threshold provides some protection against over-subtracted skies in the mea- 
surement of the disk scale length. Similar tests were performed with our 2D decomposition 
algorithm which confirm, once again, the results above. 



A number of recent studies have described the variation of bulge shapes as a function 
of Hubble type (Andredakis et al. 1995; Moriondo, Giovanardi, & Hunt 1998; Khosroshahi 
et al. 2000; MoUenhoff & Heidt 2001), going so far as suggesting that precise values of n (i.e. 
±0.1) could be determined (Graham 2001). To our knowledge, no study to date has tested 
the reliability of the recovery of the Sersic n parameter. In order to test the sensitivity of 
the decomposition to the full range of bulge profile shapes we use mock luminosity profiles 
with values of the Sersic n parameter ranging from n = 0.2 to n = 4.0. The suite of profiles 
used all combinations of rg = 0.8, 1.5, and 2'.'5, /Xg = 18, 20, and 22 mag arcsec^^ and seeing 
FWHMs of 1.5, 2.0, and 2'.'5. The profile fits used initial estimates of n = 0.4, 1, 2, and 
4, and correct initial estimates for re, //g, and the disk parameters. The seeing FWHM was 
fixed to the correct model value. The results for the n — 1 initial estimate are presented in 
Fig. 6, where we plot the average relative fit error on n (for the 100 profiles with the same 
simulated parameters) where 



3.3.6. Sersic n Tests 




(22) 



'^model 
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versus the model n (the dashed hne at An = indicates a perfect recovery of the model 
n parameter.) Each panel shows one particular combination of /ig and r^, and the panels 
are arranged snch that the B/D ratio, for a given value of n, decreases from top to bottom 
and right to left. The different seeing FWHM values are represented by three point types: 
circles, triangles, and squares for seeing values of I'.'S, 2'.'0, and 2'.'5 respectively. 

Fig. 6 reinforces that the bulges of even nearby late-type spirals are small and not 
sampled at high enough spatial resolution to yield a stable, robust solution for n as a floating 
parameter. Given the correct value of n as an initial estimate (along with the correct initial 
estimates for the other four parameters), the algorithm normally finds the correct value 
of the model n, but any departure from the model value even by a small amount, yields 
significantly different solutions for n, or the fit may simply fail (as indicated by the vertical 
lines in the figures). For most of the parameter combinations, an offset of ~50% in the initial 
estimate of n yields a ~50% error on its determined value. 

Similar tests using the 2D algorithm show a shghtly more robust recovery of the model 
n parameter based on incorrect initial estimates, but the recovery efficiency is still poor and 
results based on a fioating initial estimate of n are questionable. We are here faced with an 
under-determined optimization with too few independent data points for too many model 
parameters (at least three for the bulge). The strong covariances between n, He, and rg 
(<7n,^e) '^n,re) prcvcut a uuiquc determination of n with this NLLS code. We actually find the 
best-fit n by grid search, holding n as a fixed parameter, solving for a range of values, and 
using the xfn defined in § 3.1 to determine the best fit. Further simulations showed this 
technique to be fully reliable for the bulges considered here. 

It is difficult to estimate the error on n. Either we use a grid search and n is fixed, or 
n is kept as a fioating parameter and varies widely given wrong initial estimates. Based on 
a test with fioating n but correct initial estimates for bulge parameters, typical seeing and 
sky errors modify n by no more than 20%. 

Based on 2D B/D decompositions of simulated images and JHK images of 40 bright 
spirals, Mollenhoff & Heidt (2001) estimate that recovery errors for all the standard fit 
parameters (J^, /i, /g, Tg, and n) are less than 15%, comparable to our findings. They also 
tested for variable estimates of the sky level and seeing width, though no clear description 
of their technique is presented. 
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3.4. Summciry of the Simulations 



The tests performed in § 3.3, for idealized galaxies, allow us to define a set of guidelines 
for the reliability and limitations of our 1D/2D decompositions: 

• Initial estimates for bulge and disk parameters are unimportant provided that 



• Seeing errors must be accounted for in all bulge parameter studies. For ^ FWHM 
and a seeing measurement uncertainty at the 35% (15%) level, the bulge Ve can be 
trusted to within 10-25% (0-10%), and i^e to within ± 0.1-0.4 (0-0.2) mag arcsec'^ . 
For Te ~ FWHM+ 1 the errors improve to within 0-15% (0-10%) for r^., and ± 0.0-0.2 
(0.0-0.05) mag arcscc^^ for fig. There is no appreciable effect due to seeing on the disk 
parameters (less than 1%). 

• Sky subtraction errors dominate disk parameter errors (~5-15%) and are non-negligible 

(up to 25%) for bulges whose effective surface brightnesses are less than one magnitude 
brighter than the sky brightness (i.e. for {/isky — A'e} ^ 1)- 

• The sampling of late-type nearby bulges may not be high enough to constrain the 
Sersic n exponent uniquely as a free parameter. Iterative model fitting schemes should 
be tested for this. Our approach uses a grid search. 

• Typical seeing and sky errors modify n by no more than 20%. 

• The 2D decomposition technique does not provide a significant improvement over the 
ID method for the recovery of axisymmetric structural parameters to warrant the extra 
computational effort. 



re > (0.3)^/" * FWHM 



and 




-0.75n+1.15 forn<1.0 
0.2 n + 0.2 forn>1.0 



Armed with these basic guidelines we can now turn our attention to real data decom- 
positions using the ID technique. 
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4. Bulge-to-Disk Decompositions 
4.1. Outline 

Our study of structural properties and the variation of galaxian parameters as a function 
of wavelength uses the multi-band (BVRH) data set of late-type spiral galaxies of Courteau, 
Holtzman, & MacArthur (§ 2; Paper II). Most galaxies have at least one set of BVRH images, 
and we use multiple observations for 54 galaxies to estimate systematic errors. Other B/D 
decomposition analyses have used larger samples (e.g. Baggett et al. (1998)) but lack the 
crucial multi-wavelength information. 

We aim to develop a stable and versatile prescription to characterize structural evolution 
of the bulges and disks of galaxies. However, just as any morphological description of galax- 
ies (e.g. Hubble types) depends on the wave band, intrinsic structural parameters are also 
expected to vary with wavelength due to stellar population and dust extinction effects. Thus, 
multi-wavelength information is required for any accurate description of galaxian structural 
parameters. 

Physical differences in the shape and size of bulges among galaxies are also expected 
depending on how they were formed. Formation by accretion processes (e.g. major/minor 
mergers) can account for steeply rising de Vaucouleurs light profiles in the central parts of 
galaxies (e.g. van Albada 1982), while secular evolution would yield exponential distribu- 
tions, with or without a core, of the central light. The formation of small bulges is indeed 
largely attributed to secular processes and redistribution of disk material (see § 5.3). The 
present study is a natural extension of de Jong's (1996a) structural analysis of 86 face-on 
spirals with BVRIHK imaging, and Graham's (2001) re-investigation of de Jong's data. 

De Jong's ID and 2D B/D decompositions estabhshed significant parametric variations 
at different wavelengths. Given the intrinsic limitations of the data modeling (i.e. over- 
determination of the parameter space), his B/D fits also used a fixed Sersic n parameter (see 
§ 3.3.6), but limited to values of n = 1, 2, and 4 bulges. De Jong's analysis, and that of 
Courteau et al. (1996) who performed ID profile decompositions for 290 r-band luminosity 
profiles, supported the notion of exponential bulges and disks and a tight correlation of 
B/D scale parameters in late-type spirals. Evidence for this correlation was challenged by 
Graham & Prieto (1999) but later validated by Graham (2001) who re-modeled de Jong's 
thesis sample with a ID B/D decomposition technique^. His results support a range in the 
Sersic shape parameter from large (n ~ 2 — 3) to small (n > 0.5) values for early- to late- 



Graham's B/D analysis uses an unconstrained (floating) Sersic shape parameter (see § 4.5.1). 
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type spirals. Aware of the inadequacy of basic B/D decompositions in fitting the bulge shape 
parameter due to poor data resolution and strong covariances with other bulge parameters 
(de Jong 1996b; Broeils & Courteau 1997, see § 3.3.6), we were compeled to revisit this 
issue with our own well-tested technique and a more extensive data base. 

Our approach involves B/D decompositions with fixed n values that sample the full 
parameter space of spiral bulges, from n = 0.1, 0.2, 4.0. The fit solutions are filtered 

out on basis of relative ^-nd criteria based on our simulations (§ 3.4). We are concerned 
below with the derivation of robust B/D parameters for each galaxy profile. We compare 
our results with Graham (2001) and others, and test for any B/D parameter correlations in 
§ 5. 

The following is based exclusively on results from ID B/D decompositions. The fact 
that we do not model non-axisymmetric shapes (bars, rings, oval distortions) lessens the need 
for more computationally intensive 2D B/D decompositions, as our simulations showed no 
improvements using the 2D over the ID decomposition method for axisymmetric stucture. 



4.2. B/D Initial Estimates 

In order to determine the range of best fitted bulge and disk parameters, we need to assist 
the minimization program in finding the lowest possible (data— model) x^- From analysis of 
our mock images and profiles, we have found that any reasonable initial estimates for the disk 
parameters yields a robust solution. We base our initial estimates for the disk parameters h 
and /Iq on the "marking the disk" technique, where the linear portion of a luminosity profile 
is "marked" and the selected range is fit using standard least squares techniques to determine 
its slope. Clearly, the resulting fits are very sensitive to the adopted baseline. We tested 
various choices for the fit start and end points for our galaxy profiles including: full profile 
fit, starting points of 0.2r„iax and OAr„iax out to r^ax, and a fixed baseline shifted along 
the length of the profile and tracking 8 different locations. Additionally, we also tested the 
"moments method" of Willick (1999). The discrepancies between the different fits are large; 
(~10% on average and up to ~100% for the worst cases), but we found that the 0.2rmax 
to Tmax bascline yielded the most reliable fits (as judged by eye). The inner boundary is 
chosen to exclude the major contribution of a putative bulge or Type-II dip and Vmax is the 
radius at which the surface brightness error has systematically reached values greater than 
0.12 mag arcsec"^ (beyond which the data become too noisy to be trusted). The fits using 
the 0.2rmax to rmax baseline provided fits that were more than adequate as initial estimates 
for the disk parameters in the decompositions. 
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Flexibility in the choice of bulge initial parameters is, however, only afforded outside 
a certain range of bulge sizes relative to the seeing disk. Moreover, observed galaxy pro- 
files show significantly more variety than the idealized profiles from which these conclusions 
were drawn (e.g. we did not model Type-11 galaxies, or the presence of strong spiral fea- 
tures). Accordingly, we explore three different sets of initial bulge parameter estimates to 
protect against local minima in the parameter space. Initial bulge effective parameters were 
determined from: 

• Subtraction of the disk fit (based on the "marking the disk" technique) from the original 
profile leaving only the bulge light, rg is then computed non-parametrically from the 
data by summing up the light up to the radius which encloses half the total light of 
the bulge. Thus //e = l^ife)- 

• Tg = 0.15/1 and lie = iiq, where h and /iq are determined from the "marking the disk" 
technique. 

• (6„/log(e)) * 0.15/i and //e = {bn - bn=i) + A«o- 

The second set of initial estimates was motivated by the bulge/disk structural correlation 
found by Courteau et al. (1996) Wc added the third set of initial estimates which attempt 
to scale Te and fi^. more appropriately to the different values of n. No specific set of initial 
estimates worked better for all cases, though the S*""^ method may be the least attractive. 
It failed to provide reliable solutions (i.e. the fit failed or the values were large) in most 
cases, but in a few cases it also yielded the only viable solution. 

4-2.1. Seeing and Sky Treatment 

"Bulges" of late-type spirals arc small and their luminosity profiles can be severely 
affected by atmospheric blur. In principle, if the blurring from the atmosphere (seeing) can be 
measured accurately, it can also be corrected by Fourier deconvolution. In practice, however, 
deconvolution amplifies noise, and the seeing FWHM is subject to measurement errors. In 
§ 3 we used extensive simulations with a wide range of input parameters and various values of 
n to derive a space of recoverable parameters under specific seeing conditions, accounting for 
the typical measurement errors of our data. Seeing is accounted for by convolving the model 
profiles with a Gaussian PSF (Eq. 6) whose dispersion is measured from field stars. In order 
to account for seeing measurement errors, each profile is modeled with three different values 
of the seeing FWHM: the nominal measured value and ±15% of that value. A mean seeing 
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uncertainty of ±15% was used rather than the individual errors per measurement as these 
fluctuate greatly due, in large part, to the different number of stars in each measurement. 

Sky subtraction errors, of order ~ 0.5-1.0% in the optical and ~ 0.02% in the H-band, 
were also examined carefully (§ 3.3.5). The sky brightness measurement error is accounted 
for in B/D decompositions by using three different sky levels: the nominal measured value 
and ±0.5% (optical) or ±0.01% (H-band), of that value. 

Each profile is thus reduced 3 times for each different combination of and //g initial 
estimates, times 3 seeing FWHM values, times 3 sky values, and times 40 different fixed 
values of n, for a total of 1080 decompositions per profile. 



The 1080 decompositions for each profile are first vetted on the basis of structural 
criteria determined from our simulations (§ 3.3). A decomposition is deemed acceptable if 
it meets the following criteria: 



• B/D < 5 

• h < 15 kpc ; Te < 50 kpc 

• Te/h < 1 

The first constraint is derived from our simulations and effectively eliminates small 
bulges whose sizes are comparable to, or smaller than, the seeing disk. The remaining 
constraints are based on physical considerations and help eliminate solutions with small 
values but unrealistic parameters for late-type galaxies. Note, however, that these physical 
constraints are rather generous and do not contribute any subjective bias. 

The successful decompositions are then ranked on the basis of two indicators: (a) a 

global x^, X^h computed for the full SB profile from r = to Vmax] and (b) an inner x^, 
Xin, which includes only the central regions of the galaxy from r = to twice the radius, 
rb=d, where the intensities of the fitted bulge and disk are equal (see § 3.2). Xin adopted 



4.3. Data Filtering 



• re > (0.3)^/" * FWHM and 




for n < 1.0 
for n > 1.0 
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to increase the sensitivity of the goodness-of-fit indicator to the bulge area^. The radius 
rf,=rf is clearly a function of the bulge shape and may change from small to large n (see e.g. 
Fig. 2). Thus we use a per degree of freedom to remove any dependence of the normal 
to a changing r^. Because of the presence of spiral arms and other non-axisymmetric 
features which we do not attempt to model, the reduced is always, in principle, greater 
than unity. However, some of our solutions may have values less than unity indicative of 
an over-determined system (correlated parameters), or over-estimated errors. 

We first rank the solutions according to their and preserve only the better half. The 
reduced set is then ranked according to values and the bottom half of the distribution 
is discarded. This process is iterated at least twice, or until we reach 50 or fewer solutions. 
Solutions with x% greater than 50 in this final subset are discarded. 

Ideally, the minima for the distributions of x% and xin values should agree to a common 
value of n, but differences may exist. We search the final < 50 solutions for a common 
solution, starting at the minima of each x^ distribution. If the n values corresponding to 
the two minima do not agree, the n values for the next smallest x^ values are compared 
(with the lower x^ values and with each other) , and this process is iterated up to three times 
until a match is found. If this process did not converge, i.e. there is no true minimum in 
the ((Xgi) X%i) - ^ space), a final solution is chosen corresponding to the minimum value of 
(Xgi/min(Xgi) + xL/min(xL)). 

Figs. 7 & 8 show examples of the distributions of (x^^, xD versus n where XpZoftoi' = 
Xgi/min(Xg,ji,J and xlnner' = xL/min(x-„,/iit), where min(xjiit) is the minimum x^ value 
from the set of (< 50) filtered solutions. Note that these minima do not necessarily cor- 
respond to the lowest value of the respective distributions from all 1080 solutions, as the 
initial absolute minima may have been filtered out (i.e. a poor combination of x^^, xln ^'^'^ 
a given solution). Thus, the normalized x^'s may be less than one (as is easily seen in the 
leftmost plot of Fig. 8). In these plots, the left panels show the x^ distributions for all 1080 
decompositions, while the right panels display only the < 50 solutions remaining after the 
iterative filtering scheme described above. 

Fig. 7 highlights the sensitivity of our technique for two V-band observations of UGC 929 
taken under different seeing/sky conditions. The left figures show a fairly well-behaved 
solution favoring n = 0.6 and the figures on the right plot show a rather messy solution 
favoring n — 0.8. The seeing conditions were worse and the sky was much brighter for the 
observation shown on the right which could explain the noisy distributions of both the x^; 



^Notc that our algorithm minimizes the x^; only. The Xj„ is calculated and used as a discriminator only 
after the algorithm has converged. 
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Fig. 8 shows two different behaviors of xli ^oi profiles with very well-behaved Xin- The 
plot on the left for our UGC 784 B-band profile illustrates the need for an additional, 
more discriminating statistic for the bulge region. Decompositions based solely on the Xgi 
goodness-of-fit indicator may result in fits, hke the one shown on the right side of Fig. 9 
(dashed-dotted bhic line). However, Fig. 8 a) clearly shows that the fit with n — 0.6 is a far 
superior match to the bulge shape, as indicated by the Xin behaviour. 

The final step of our filtering procedure entails a visual inspection of the final decom- 
positions. The criteria for user examination include information from multiple exposures 
and multi-band reductions for a given galaxy. Profiles and/or solutions with the following 
pathologies were eliminated from the final sample: 

• disk profiles that are too short for proper fitting 

• no obvious, extended, underlying exponential structure for the disk (occurs predomi- 
nantly in Type-II profiles) 

• unphysically large fitted bulge 

• unrealistic disk fit for Type-II profiles. The fit is tipped below the true disk to account 
for the Type-II dip near the bulge-disk transition region leading to scale lengths that 
are biased high (e.g. see Fig. 12 for UGC 12527 for examples of "bad" fits which were 
eliminated from the final sample) 

• large deviations between solutions for multiple observations of a given galaxy. 

Not surprisingly, most of the ehminated profiles are Type-II systems. We caution that 
even the Type-II profile decompositions that survived the full sorting process may not provide 
the ideal description of their complex surface brightness distributions. Cleary these Type-II 
profiles cannot be properly modeled with just a Scrsic bulge and exponential disk. Out of 
523 images/profiles, a total of 341 passed our acceptance criteria. 

4.3.1. Preferred Sky and Seeing 

Histograms of the preferred seeing FWHM and sky offsets for all decompositions in all 
four bands are shown in Fig. 10. Typically, a lower sky brightness level is preferred by our 
algorithm. In some cases, this could be explained by an over-estimated sky level, but it may 
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also be due in part to profiles with truncated outer disks as in Fig. 13. Our program prefers a 
slightly larger seeing FWHM than measured. This could be the result of an under-estimated 
FWHM, or perhaps differences between the idealized Gaussian model and the real seeing 
PSF. Solutions with variable sky/seeing estimates were retained in the final solution set for 
assessment of parameter errors. 

4.3.2. Effect of 

I max 

Of significance to the fit results is the maximum radius used in the decompositions. We 
have used the full observed profile out to radii where the surface brightness errors reached 
above 0.12 mag arcsec^^ . There is no absolute definition to the edge, r^ax^ of a disk 
and a different selection could yield different results. To test the sensitivity of our parameter 
determinations to the chosen value of Tmaxi we rc-dccomposed the profiles as described above, 
but with a fit baseline extending only to 0.75 x r„iax- A comparison of the results from the 
two techniques, prior to eyeball filtering, shows good agreement and we chose to keep the 
larger basehne to avoid discarding good data. 

4.4. Decomposition Examples 

There is no room for a full display of our catalog of final decompositions, but a few 
examples are shown in Figs. 11-14. The full catalog of decomposition plots is available upon 
request from the authors. 

In these figures, the solid black circles are the data points, the black dots show the 
sky error envelope (from the measured sky error), the dashed and dashed-dotted lines show 
the disk and bulge fits respectively, and the sohd line is the total (bulge+disk) fit. The fits 
are all seeing-convolved using the best selected seeing values. The bottom panel shows the 
fit residuals where A//(r) represents the (data— model) . Fig. 11 shows an example of the 
quintessential Type-I profile at all wavelengths. Fig. 12 shows a Type-II/Transition galaxy 
whose Type-II signature significantly weakens from the optical to the infrared. Fig. 13 shows 
a Type-I profile with an outer truncated disk. Such decompositions will presumably favor 
an under-subtracted sky in attempt to align the inner and outer parts of the disk. Here is 
an example where our procedure with an over/under-estimation of the sky and an infinite 
exponential disk model may not be adequate since the outer disk truncation appears real (as 
detected in all four bands). We also show an example of a nearly bulgeless system in Fig. 14. 
Our sample is divided into 52 Type-I, 53 Type-II and 16 transition systems, of which 18 
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truncated and 7 bulgeless disks are identified^°. 

4.5. Distribution of the Sersic n parameter 

Fig. 15 shows histograms of the Sersic n parameter for all the final fits (left) and good 
fits only (right) after user-examination as described above. The distribution of n has a 
definite range, implying that not all late-type bulges are best described by an exponential 
profile, but the mean value is very close to one. This result agrees with Graham (2001) and 
recent N-body simulations of galaxy evolution (§ 5.3). 

4.5.1. Floating Sersic n 

In § 3.3.6 we showed that resolution limitations prevented stable fitting of the Sersic 
shape parameter n as a free parameter. To illustrate the effect a fioating n can have on fitted 
parameters we re-decomposed all of our galaxy profiles leaving n as a free parameter (e.g. 
akin to Graham 2001). The results are shown in Fig. 16 for three different initial guesses for 
n (0.2, 1.0, and 4.0). 

The histograms of the resulting distributions of n reveal a strong bias towards the 
adopted initial estimate. All 3 distributions show a large peak at n = 0.1, indicative of poor 
bulge fits. The histogram for the n = 1.0 initial estimate looks somewhat similar to our own 
constrained solution (Fig. 15), but this is somewhat fortuitous given the closely-exponential 
nature of spiral bulges. Note also the non-Gaussian tail in Fig. 15 is not reproduced in 
Fig. 16 for the n = 1 initial estimate case. Fig. 17 shows a comparison of the values from 
the fioated versus fixed n solutions. (When no suitable fit was found, all parameters were set 
to as indicated by the points lying on the axes; note the large number of fit failures in the 
fioated n case.) Note the large discrepancies in Xin X^i between the two methods. Thus, 
while the final distributions for the n — 1 initial estimate and our constrained n procedure 
look similar, significant differences may exist between individual decompositions. 



^°The fraction of a given galaxy class should be interpreted with care since our sample is not volume- 
limited. 
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4.6. Error of a Single Measurement 



An important feature of any decomposition technique is the stabiUty of the final results 
for repeat observations of a given system. Our sample has 50 profiles for which multiple (two 
to four) observations exist, allowing for a direct measure of the reliability of our decompo- 
sitions. Table 1 gives the mean and mean standard deviation of the five model parameters 
from repeat observations with, 



where x is the fit parameter, n is the number of observations for a given profile, and N is the 
number of profiles with repeat observations. The average errors from repeat observations of 
Type-I profiles are ± 14% for n, ± 0.2 mag arcsec^^ for fie, ± 13% for r^, ± 0.05 mag arcsec^^ 
for //-o, and ±3% for h. Clearly, determinations of disk parameters are much more stable 
than for bulges. Error terms quoted below correspond to the l-cr deviation, unless otherwise 
noted. 



The overlap between our sample and de Jong's thesis catalog (de Jong & van der Kruit 
1994) amounts to only 3 galaxies. Direct comparison of our SB profiles shows excellent 
zero-point and overall shape agreement (Paper II); however, our B/D decompositions differ 
somewhat, as shown in Table 2 (note that de Jong uses fixed n — 1). Also shown in that table 
are decomposition parameters for the same galaxies by Graham (2001) (same data as de Jong, 
hut using a range of n). We find scale length differences at the 10% level with de Jong and 
Graham, consistent with, or slightly better than, typical variations between different authors 
(Knapen & van der Kruit 1991). A comparable dispersion is measured between the scale 
lengths of Graham and de Jong based on 82 R-band profile decompositions. Graham's scale 
lengths are, on average, smaller for small galaxies and larger for big galaxies (apparent size) 
than de Jong's. We find systematically larger disk scale lengths than de Jong (based on only 
7 comparisons.) We verified that sky under/over-estimates cannot account for any difference 
with de Jong. De Jong's algorithm gives more weight to the outer part of the disk, possibly 
explaining the shorter disk scale lengths. For profiles with outer truncated disks or Type-II 
decrements, greater weight in the outer parts favors the outer disk curvature and thus steeper 
disk fits. 

Bulge parameters between us and Graham match reasonably well for the first two galax- 
ies but differ substantially for UGC 3140. It is however difficult to establish trends based 
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4.7. Compcirison with Other Authors 
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on just 3 comparisons. We can, instead, broadly compare our respective distributions of 
Sersic n with morphological type. This is done in Fig. 18 for comparison with Fig. 10 of 
Graham (2001). The general features are similar, but we find a wider range of n values for 
the later-types possibly due to the larger number of Scd/Sd galaxies in our sample. Another 
favorable comparison of bulge parameters with Graham is shown in Fig. 20 (see § 5.2.1). 

We also find an overlap of two galaxies, NGC 3512 and NGC 7782, with the sample 
of Baggett et al. (1998). As with de Jong and Graham, disk parameters agree within 10%. 
Bulge parameters from Baggett et al. are missing for NGC 3512, and those listed for NGC 
7782 (both V band) differ quite substantially from ours. These authors find jif. = 10.88 and 
Te = 0.2 for a de Vaucouleurs bulge and wc have /Xe = 20.1 and Ve = 3.2 for n=l (best 
fit) or He = 25.2 and Ve — 83.8 for n—A (very bad fit with high reduced x^)- Note that 
seeing estimates were comparable. Surprisingly, their /i^ is nearly 10 magnitudes brighter 
than their //q = 20.3 (we also find //q = 20.3)! We find this pathology in a number of their 
bulge decompositions (see e.g. their Fig. 2) where the models often overshoot the data at 
the center. 

We conclude this section by noting that disk scale lengths between us and other authors 
differ at the 10% level. Our bulge parameters are also qualitatively consistent with those of 
Graham. 



5. Discussion 

Simulations of galaxy profiles and images (§ 3) and careful B/D decompositions (§ 4) 
have led to a final set of structural parameters for late-type spiral galaxies (Table 4 in Ap- 
pendix B) . These data can now be examined for intrinsic structural variations and sensitivity 
to dust and stellar population effects. The outline of this section is as follows: First, we ver- 
ify in § 5.1 that our solutions are not affected by projection effects. Wc then discuss in § 5.2 
B/D parameter variations both in the context of profile type differences and wavelength de- 
pendence. In light of existing limitations in our modeling of Type-11 profiles, our conclusions 
will be based mostly on properties derived from Type-I profiles. These will enable us to 
examine the viability of secular evolution models for disk galaxies (see § 5.3). 

5.1. Inclination Dependence 

In order to test for projection effects, we plot the distributions of fi^ and r^, as well as 
disk iiQ and /t as a function of ellipticity, s = 1 — b/a, in Fig. 19. The surface brightnesses 
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are only corrected for Galactic extinction and cosmological dimming (as in § 2); thus the 
/xo and fie values should be considered as upper limits (i.e. effective brightnesses are too 
low). No trends with ellipticity are seen, including the Sersic n parameter and ratio of disk 
scale lengths (not plotted). Furthermore, Types 1, II and Transition arc not confined to any 
particular inclination range showing that the Typc-11 phenomenon is not an accentuated 
feature due to line-of-sight extinction (e.g. Type-II galaxies are not preferentially inclined 
with the plane of the sky) . 

5.2. Bulge/Disk Pcirameters 

Tabic 3 shows the range of fitted parameters at BVRH wavelengths for all galaxy profile 
types (Typc-1, 11, and Transition). The number of Type-11 and Transition systems included 
in this table (e.g. only 4 decompositions for Transition galaxies in the B-band) is drastically 
reduced from our original distribution as many of them did not pass our validity criteria 
(§ 4.3). Note that the parameters for Transition profiles at H-band broadly match those of 
Type-I's at that wavelength. 

The Sersic shape parameter n for Type-I galaxies is near unity, within the errors, for 
all wavelengths. Thus, we advocate that the natural, intrinsic distribution of the Sersic 
n parameter for late-type spirals has a mean near 1.0 (with (T„ ~ 0.4; see Fig. 15). By all 
accounts, bulges of late- type spirals are well- approximated, on average, by a pure exponential 
(luminosity /mass) density distribution. 

The distributions of Sersic //g and are broad, indicative of the range of bulge types in 
our sample. Effective radii are typically less than 1 kpc. Those of Type-II profiles are even 
smaller and seemingly better determined than Type-I's but this is predominently an artifact 
of our limited 2-component modeling. Examination of Type-II profile fits shows that the 
model disk is typically shallower (than the "true" disk), as the fit accounts for the fainter 
bulge/disk transition dip, and bulge effective radii are naturally confined to a smaller range. 

The distributions of disk scale lengths and their ratios show a clear decreasing trend 
as a function of wavelength (as noted by de Jong (1996b)). This statistically significant 
effect, detected for all profile types, can be explained either by a high concentration of older 
stars and/or dust in the central regions of the disk. Absorption by dust alone can account 
for the scale length ratios that we measure (see e.g. Evans (1994), Fig. 5). Evans' models 
do not consider scattering but for the nearly face-on galaxies considered here, its effects 
are negligible (Byun et al. 1994; de Jong (1996c)). The color gradient analysis of de Jong 
(1996c) using stellar population and dust extinction models suggests, however, that dust and 
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metallicity play a minor role but that age is be the dominant factor. Preliminary analysis 
of our photometric data with the latest stellar evolutionary and dust models (MacArthur et 
al., in prep.; hereafter Paper III) suggests a combination of effects. The interpretation of 
color gradients is non-trivial and may ultimately require a full spectroscopic investigation to 
convincingly disentangle the effects of age, dust, and metallicity. 

5.2.1. B/D Scale Ratios 

In Fig. 20, wc plot rg vs. h for our Type-I decompositions (solid symbols) and those of 
Graham (2001) for dc Jong's BRK Type-I SB profiles (open symbols). This figure provides 
the basis for a renewed discussion of the suggestion by Courteau ct al. (1996) of structural 
coupling between the bulge and disk of late-type galaxies. The large dispersions in the 
Tg and /i'^ (see Table 3) nearly cancel out to yield tighter Ve/h correlations. For Type-I 
profiles we find {ve/h) ~ 0.22 ± 0.09 at all wavelengths, corresponding to (/ibuige/^disk) = 
0.13 ± 0.06 for n — 1. This result is also borne out in the H-band Transition profiles (see 
Table 3). For comparison, Courteau et al. (1996) found (/ibuige/^disk) ~ 0.10 ± 0.05 (or 
{ve/h) = 0.15 ± O.OS)-*^^. The latter result is also in agreement with Graham (2001) who 
finds {Ve/h) = 0.2 (no quoted dispersion, but it is somewhat larger than ours judging from 
Figs. 20 and 21) for early and late-type spirals in the K-band. This is consistent with a 
scenario where bulges of late-type spiral galaxies are more deeply embedded in their host 
disk, than earlier-type bulges. In such an "iceberg" scenario (e.g. Graham 2001), bulges and 
disks can preserve a nearly constant Ve/h but show a great range of He for any given Ve- 
in Fig. 21, we show Ve/h as a function of morphological type from our (solid symbols) 
and Graham's (open symbols) decompositions. A mild trend with Hubble type is seen with 
(rjh) = 0.20 - 0.013(T - 5) {la = 0.09), ranging from (re/h) ~ 0.20 for late-type spirals to 
(re/h) ~ 0.24 for earher types. Comparison of our and Graham's decomposition parameters 
in Table 2 shows that large deviations may exist, thus only our data points were included 
in the fit of (rg/ h) vs T above. Data for different bands scatter evenly about the mean hne. 
More data at earher and later types would be needed to firm up this trend. It is nonetheless 
remarkable that early and late-type systems are described by very similar scaling relations, 
thus suggesting comparable formation and/or evolution scenarios. 



^^Thc study of Courteau ct al. (1996) combined the r-band decompositions of Broeils & Courteau (1997) 
and the K-band decompositions de Jong's (1996a) thesis study but no distinctions were made between Type-I 
and Type-II profiles. 
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5.3. Test of Seculeir Evolution in Late- Type Spirals 

This work has provided confirmation of two important structural signatures of spiral 
galaxies which must be addressed by models of structure formation: 

• The underlying surface brightness distribution of late-type spirals has a range for the 
Sersicn parameter from 0.1-2, but is best described, on average, by a double- exponential 
model of bulge and disk, such as found in Type-I profile galaxies. 

• Bulges and disks of late-type spirals are coupled, with (re/h) — 0.22±0.09; or (/ibuige/^disk) 
= 0.13 ± 0.06; at all wavelengths. A mild trend with Hubble type is also detected with 

a range (re/h) ~ 0.20-0.24, from late to early-type spirals. 

The first result describes the large-scale appearance of bulges. Analyses of HST images 
have shown that a significant fraction of bulge nuclei have power-law profiles (r < 500 pc; 
e.g. PhiUips et al. 1996; Balcells 2001) and host a central compact source (CaroUo 1999). 
The extent of these nuclear sources (< 0'/3 for cz < 2500 kms~^) is smaller than our images' 
pixel size and smoothed out by seeing. We thus ignore their effects on the bulge fight profile 
in this analysis, but caution that our bulge parameters are to be considered upper limits if 
a significant nuclear component is present. 

A natural interpretation of the near constancy of B/D size ratios in late- type spirals is 
that their bulges formed via secular evolution of the disk. This scenario is possible if disks are 
bar-unstable, which can be triggered by the global dynamical instability of a rotationally sup- 
ported disk or induced by interactions with a satellite and if significant angular momentum 
transport is feasible (e.g. Martinet 1995; Combes 2000; see the collection of papers in Carollo, 
Ferguson & Wyse 1999 for comprehensive reviews). For bar-unstable disks, in particular to 
vertical deformations, the inner disk material is heated up to 1-2 kpc above the plane into 
a "bulge" via resonant scattering of the stellar orbits by the bar-forming instabifity. This 
in turn, catalyzes funneling of disk material into the central regions and generates outward 
transport of disk material in the outer parts. Gas fiows must also be invoked to explain the 
higher spatial densities of bulges compared to the inner disk. Such a model is expected to 
produce correlated scale lengths and colors between the disk and its central regions, as ob- 
served (e.g. Terndrup et al. 1994; Peletier & Balcells 1996; Courteau 1996b). A "bulge-like" 
component with a nearly exponential profile is expected from non-axisymmctric disturbances 
that induce inward radial flow of disk material (Pfenniger & Friedli 1991; Zhang & Wyse 
2000, and references therein). The longer the disk-bar heating interaction, the greater the 
extent of the disk exponential profile (Valenzuela & Klypin 2002). The evolving exponential 
attractor is an empirical result well established in simulations, but it lacks, at present, a 
theoretical explanation (Pfenniger 1999). 
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Although a bar can grow spontaneously ( ^ 20 Myr) from small scale fluctuations in the 
inner disk, an external finite perturbation can catalyse its growth. However, coUisionless 
mergers seem unsuited to growing the exponential bulges of present-day late-type spirals, 
though they may contribute to the increase in Scrsic n parameter seen toward earlier types in 
proportion to the accreted satellite mass (Barnes 1988, Aguerri et al. 2001). The spontaneous 
or triggered formation of bars also suggests that the Hubble type of galaxies can change well 
after the formation of the disk (Pfenniger 1999). All of our bulges are smaller than a disk 
scale length and could be created by purely bar-related processes. Instead, accretion of 
galaxy satellites is required to make bigger bulges, either before or after formation of the 
host disk. 

Secular evolution models of stellar and gaseous disks, especially through cosmologically- 
motivated three-dimensional N-body simulations, have seen significant developments in the 
last decade. For example, the cold dark matter (CDM) hierarchical hydrodynamical simula- 
tions by Saiz et al. (2001) and Scannapieco & Tissera (in prep.) show that secular processes 
can occur naturally during the formation of spiral disks and play an important role in the 
regulation of star formation and the determination of the dynamical and structural proper- 
ties of these systems. On average, the simulated disk systems are shown to be characterized 
by a double exponential profile which naturally emerges within the hierarchical clustering 
scenario. These results are based on a stellar formation process implemented in such a way 
that it succeeds in forming compact bulges that stabilize disk-like structure allowing the con- 
servation of an important fraction of their angular momentum during the violent phases of 
their assembly. Fig. 22 shows the distribution of final Sersic n parameters for relaxed present- 
day late-type disks by Scannapieco & Tissera (in prep.); their models reproduce our results 
(Fig. 15) very nicely. The double-exponential structure of bulge and disk may not always be 
the final relaxed state of an object, but whenever n ~ 1, the B/D scale ratio (/ibuigc/^disk) 
takes its nominal value of 0.15. These models do not have bulges with n < 0.7, possibly due 
to limited resolution and/or excessive angular momentum transfer that supernova feedback 
could help prevent. 

Simulations by Pfenniger (2002; private comm.) of self-gravitating disks forming bars 
which may later dissolve into a bulge-like component also show a nearly universal ratio Tg/ h, 
in agreement with observed values, which is related to the stellar dynamics of the barred 
system (i.e. relative position of the vertical to horizontal resonances). The bar length is 
related to the initial rising part of the rotation curve (yielding a scale) , and the corotation of 
bars is proportional to their length. The corotation fixes the positions of the other resonances, 
which in turn fix the maximum extension of bulges made from resonant heating, as indeed 
the vertical resonances are strong only within the bar. This mechanism would thus set the 
upper limit for the allowed range in rg/h. 
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The A'"-body simulations of Aguerri et al. (2001), which consider the growth of galactic 
bulges by mergers, also suggest that the final B/D scale ratio {ve/h) does not scale with 
the B/D luminosity ratio. These authors show that the disk scale length h can increase 
from 15% (low mass retrograde satellite) to 65% (high mass direct satellite) while (rg/Zi) 
would decrease from 0.21 to 0.14, or 33%, in the most extreme case. One can thus infer 
{ve/h) = 0.17 ± 0.03, independent of B/D luminosity ratio, in good agreement with our 
findings. Their simulations are however limited to a small range of initial re and a more 
complete investigation with a broad range of and h values is needed to establish the 
fundamental nature of the B/D scale ratio. 



5.4. Type-II Profiles 

The above scenarios for secular evolution naturally produce the double-exponential char- 
acter of the bulge and disk radial luminosity profiles for late-type systems. However, over 
half our sample of 121 late-type spiral galaxies show strong deviations from this simple 
two-component description. Other authors (Kormendy 1977; Baggett et al. 1998) have con- 
sidered inner disk truncation (plus de Vaucouleurs bulges) as an alternative to modeling 
Type-II light profiles, with 

Idiskir) = lo exp {- [r/vo + {rhoie/rT]} (24) 

where rhoie is the truncation radius and n ~ 3. As discussed in § 3.2, we do not consider this 
approach at the present, but its potential merits should not be overlooked. 

N-body simulations (e.g. Norman, Sellwood, & Hasan 1996; Valenzuela & Klypin 2002) 
reproduce Type-II surface density profiles as a result of the redistribution of central stars into 
a ring by a bar-like perturbation. Approaching the centers where the component called bulge 
and the component called exponential disk overlap, one cannot tell, in these simulations, 
if a star or particle belongs to which component. Galaxy centers may recurrently move 
from a barred to an unbarred phase and undergo continuing bulge building as the bars 
dissolve^^ (Norman et al. 1996). Thus, the paucity of barred galaxies in our sample does 
not preclude bar-induced effects as a possible explanation for Typc-11 profiles (e.g. Gadotti 
& dos Anjos (2001)). Pre-existing bars may simply have dissolved. For example, out of 8 
barred-classified galaxies in our sample, 6 have Type-II profiles thus lending some credence to 
the bar-lens scenario. On the other hand, the most strongly barred galaxies in the Shellfiow 



""^^Only progressively larger bars in the centers of exponential bulges would be allowed to form in a recurring 
scenario due to the disrupting dynamical effect of a growing nucleus (Rix 1998, as reported in CaroUo 1999). 
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sample of ~ 300 bright late-type galaxies (Courteau et al. 2000) have mostly Type-I profiles, 
indistinguishable in shape and global properties from the profiles of unbarred Type-I galaxies 
(Courteau et al., in prep.). ^From an inhomogeneous sample of 167 spiral galaxies, Baggett 
et al. (1996) find only a weak tendency for barred galaxies to have a higher occurence of 
Type-ll profiles. The link between Type-II profiles and barred galaxies is thus unsecured at 
present. 

Type-Il profiles may also be explained by extinction effects in the disk. Increased 
opacity towards the central disk can cause a depression in the luminosity profile, especially 
at shorter wavelengths. Realistic Typc-11 profiles (in shape and colors) have been produced 
with exponential distributions of stars and dust and variable layering parameters (Evans 
1994). If dust extinction causes the inner disk profile dip. Transition galaxies would just be 
a case of lesser dust content, whereas bona fide Type-II systems remain optically thick, even 
at H-band. Using far-infrared (FIR) to B-band flux ratios, and radiation transfer models for 
the dust (Gordon et al. 2001), we have tested for the origin of Type-II signature as being due 
to extinction. The FIR/B fiux ratio should be higher for the dustier systems. Unfortunately, 
our measured total FIR/B fiux ratios are statistically identical (with large scatter) for Type- 
I, Type-II, and Transition galaxies (Paper III), thus thwarting any clear interpretation. The 
IRAS 60 and 100 /xm fiuxes have too low resolution and too large errors to separate the inner 
disk dust emission from the whole galaxy. 

If stellar population effects are relevant (Prieto et al. 1992), age/metallicity gradients 
should be detected at the bulge/disk transition in Transition systems. We will further 
investigate the dust and/or stellar population origin of the Type-II dip in Paper III. 

6. Summary and Concluding Remcirks 

This study has focused on the development of rigorous B/D decomposition techniques 
using a new, comprehensive, multi-band survey of late-type spiral galaxies. We examine 
three types of SB profiles. Freeman Type-I and Type-II, and a third "Transition" class for 
galaxies whose profiles change from Type-II in the optical to Type-I in the infrared. This 
distinction is important since Type-II and Transition profiles cannot be adequately modeled 
by a simple two-component model of the bulge and disk. Thus, our main results are based 
on Type-I profiles. 

Based on extensive simulations, careful treatment of sky and seeing measurement errors, 
and repeat observations we are confident that systematic errors are < 20% for the bulge 
components, including the Sersic shape parameter, and < 5% for disk components. 
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The main conclusions from our simulations and final profile decompositions are as fol- 
lows: 

• Simulations to determine the range of acceptable solutions for any B/D decomposition 
program are crucial. The reliability of bulge model parameters is limited by the relative 
size of the bulge and seeing disk, seeing errors, the intrinsic bulge shape, sky brightness 
and errors. Disk parameters are fairly robust to systematic errors, with the exception of 
improper bulge shapes and sky errors which can have dramatic effects on both modeled 
disk and bulge components. 

• The Sersic bulge shape parameter for nearby late-type galaxies shows a range between 
n = 0.1 — 2, but, on average, their underlying surface brightness distribution is best 
described by a double-exponential model of bulge and disk. 

• Disk scale lengths decrease at longer wavelengths, indicative of a higher concentration 
of older stars and/or dust in the central regions relative to the outer disk. 

• We confirm and reinforce the result of Courteau et al. (1996) of a structural couphng 
between the bulge and disk of late-type spirals. We find (re/h) — 0.22 ± 0.09, or 
(^buige/^disk) = 0.13±0.06, independent of wavelength. A mild trend with Hubble type 
is observed with (re/h) = 0.20 - 0.013(7 - 5) {la = 0.09), ranging from (re/h) ~ 0.20 
for late-type spirals to (re/h) ~ 0.24 for earlier types. These results are consistent with 
scenarios of bulge formation in which bulges of late-type spiral galaxies are more deeply 
embedded in their host disk than earlier-type bulges. Under this "iceberg" scenario, 
bulges and disks can thus preserve a nearly constant Ve/h but show a great range of 
He for any given r^. The observed scale ratio is consistent with numerical simulations 
of self-gravitating disks and probably related to the stellar dynamics of an actual or 
pre-existing barred system. 

• The inner brightness profile signatures of Type-11 galaxies are likely explained by a 
combination of dust extinction and stellar population effects and perhaps linked to the 
occurence of a bar, but no decisive conclusion can be derived at present. 
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A. Functional form for the Sersic 6„ pcirameter 

Eq. 14 cannot be solved in exphcit closed form for 6„. Many of the numerical and ana- 
lytical solutions found in the literature agree well for n > 1 but differ significantly for smaller 
values of n. Fig. 23 shows a comparison of the two most commonly used approximations 
(short- and long-dashed curves) with the exact solution for 6„, computed to a numerical 
precision of one part in 10^ for all n < 10 (see also Fig. 2 in Graham 1999). 

As we wish to test for spiral bulges with Sersic n's as small as 0.1, we have adopted 
a formalism that is valid for all n. To maintain computational simplicity, and ensure a 
suitably accurate solution we found it practical to divide the curve into two segments. For 
all n > 0.36 we use the asymptotic expansion of Giotti & Bertin (1999) up to 0(n~^) (their 
Eq. 18), 

1 4 46 131 2194697 ^,5, 
bn^2n \ 1 H + 0(n ^) (Al) 

3 405n 25515n^ 1148175n'^ 30690717750n'i ^ ^ ^ ^ 

which is good to better than one part in ~ 10'^ in that range. However, for n < 0.36 this 
solution diverges. Due to the rapidly changing curvature in the gamma function, and thus 
bn, at small n, it would be necessary to use an unrealistic number of terms in the asymptotic 
expansion to achieve the desired accuracy. For n < 0.36, we find the best fitting polynomial 
of the form 6„ = X^I^o * where m is the order of the polynomial and the Oj are the 
coefficients of the fit given by, 

ao = 0.01945 ai = -0.8902 as = 10.95 03 = -19.67 04 = 13.43. (A2) 

This fit is accurate to better than two parts in 10'^. The wiggles in the dotted curve in 
Fig. 23 result from the polynomial nature of the fit and limited numerical precision where 
the gamma function approaches infinity. 
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B. Decomposition Results for the Type I Profiles 

Table 4 gives relevant photometric information and ID B/D decomposition results for 
the final set of Type I galaxy profiles. Decomposition results for Type-II and Transition 
galaxies are available from the authors upon request, with the caution that parameters for 
these profile types should be interpreted with care. The entries are arranged as follows: 

Column (1): (UGC number) (observation number) (passband) for each profile; 

Column (2): EUipticity, e = (l—b/a). The final ellipticity (and position angle) estimates 
correspond to an average of those values from the five contours surrounding the best isophotal 
fit in the outer disk, as determined by eye. This estimate is clearly sensitive to the presence 
of spiral arms. The typical inclination error is ~ 3 deg, independent of ellipticity; 

Column (3): Sky brightness in mag arcsec~^ , measured from 4 sky boxes located 
between the detector edges and a fair distance away from the galaxy. Typical rms sky errors, 
computed from the deviations of the mean sky counts in those sky boxes, are ~ 0.5 — 1.0% 
in the optical and 0.05% in the JR. The subscripts indicate the sky offset preferred by our 
selection process as described in § 4.2.1 and § 4.3 (and see Fig. 10), where "+" and "— " 
indicate 0.5% for optical and 0.01% for H-band over- and under-subtracted skies respectively. 
No subscript indicates that the measured sky was preferred; 

Column (4)'- Seeing FWHM values, computed as the average of the FWHMs of all non- 
saturated stars measured automatically on each image frame; typically 10 to 30 measure- 
ments per image were used for each FWHM estimate. The accuracy of the seeing estimate 
per image is roughly 20% for the optical bands and 30% for the H-band. The subscripts 
indicate the seeing offset preferred by our selection process as described in § 4.2.1 and § 4.3, 
(and see Fig. 10) where and "— " indicate 15% over- and under-estimated seeing FWHM 
respectively. No subscript indicates that the measured seeing FWHM was preferred; 

The upper and lower boundaries in the remaining columns correspond to the maximum 
and minimum values of the < 50 (out of 1080 total) solutions left after filtering (see § 4.3); 

Column (5): Best fit Sersic n bulge shape parameter; 

Column (6): Bulge effective surface brightness, yUg, in mag arcsec"^ , corrected for 
Galactic extinction and cosmological redshift dimming, as described in § 2.2; 

Column (7): Bulge effective radius, re, in arcseconds; 

Column (8): Bulge effective radius, Tg, in kpc. Converted to a physical scale using the 
Local Standard of Rest velocity, Vlg (see Paper II); 
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Column (9): Exponential disk central surface brightness, in mag arcsec ^ , corrected 
for Galactic extinction and cosmological redshift dimming as described in § 2.2; 

Column (10): Exponential disk scale length h, in arcseconds; 

Column (11): Exponential disk scale length in kpc. Converted to a physical scale 
using the Local Standard of Rest velocity, Vlg (see Paper 11); 

Column (12): Bulge-to-disk luminosity ratio, B/D, calculated using Eq. 17 in § 3.1. 
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Fig. 1. — Examples of Type- 1 (left), Type-II (middle), and "Transition" (right) SB profiles. 
Blue circles, green squares, red triangles, and purple asterisks are for B, V, R, and H-band 
respectively. The solid black lines plotted on the B-band and H-band profiles are fits to the 
outer exponential disk profile. 
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5 10 15 20 



Radius 

Fig. 2. — Sersic n profiles for different values of n. The top panel shows profiles with /ie = 21 
mag arcsec"^ and = 3'.'5 for values of n in the range 0.2 < n < 4. The table lists the 
relative light contributions of the different profiles normalized to the n = 1 case. The bottom 
panel shows the same profiles except for a constant CSB of /xq = 18 mag arcsec"^ . 
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Fig. 3. — Effect of fitting an incorrect Sersic n bulge on the disk scale length h. Each panel 
plots the average relative fitted h errors {Ah = (/;,£|.(mean) - /imodeO/^model) with solid 
symbols and connected by solid lines as a function of the model n for a bulge with =2'.'5 
and fj,e = 20 mag arcsec"^ . Red circles, green triangles, and blue squares correspond to 
seeing values of 1.5, 2.0, and 2'.'5 respectively. The three panels are for model n values of 0.2, 
1.0, and 4.0 from left to right. 
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Fig. 4. — Effect on the fitted rg value of an incorrect seeing value in the ID decomposition. 
The column plots are based on different values for the fractional seeing error used in the fit, 
where Ap-y^jjjYj = (FWHM^^gg^^ — FWHMj^^Q^gj)/FWHM^Q^g^. Each row is for a different 
value of the model FWHM, I'.'S (top) and 2'.'5 (bottom). Each panel shows the average 
relative error on Vg, Are (Eq. 21), versus the model rg. The seven curves are for different 
values of fie- 16 (dark purple), 17 (blue), 18 (red), 19 (green), 20 (magenta), 21 (cyan), and 
22 (orange) mag arcsec"^ . 
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Fig. 5. — Effect on the fitted r^. value of an incorrect seeing value in the 2D decomposition 
(compare with Fig. 4.) The solid symbols connected by solid lines indicate the average (of 
the 40 image decompositions for each parameter and initial estimate combination) relative 
error on re, Are (Eq. 21). Blue circles and red triangles are for /ie values of 18 and 22 
mag arcsec"^ respectively. Ap-yyjjjYj is as defined in Fig. 4 
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Fig. 6. — Difference between modeled and recovered values of n for a range of artificial profiles 
from n = 0.2 — 4. The Sersic exponent n is a free fit parameter and the initial estimate is 
set to n = 1. Each panel shows the average relative fitted n errors {An = (nf^^(mean) - 
^model)/''^model) ^i^h solid symbols and connected by solid lines versus the model n for 
the 9 combinations of Ve = 0.8, 1.5, 2'.'5, and fie = 18, 20, 22 mag arcsec"^ . Red circles, 
green triangles, and blue squares correspond to seeing values of 1.5, 2.0, and 2'.'5 respectively. 
The panels are ordered such that the B/D ratio, for a given n value, decreases from top to 
bottom and right to left panels. 
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Fig. 7. — Examples of xLner' (open blue triangles) and x^iobai^ (filled red squares) versus 
Sersic n distributions for the 1080 decompositions of two different V-band observations of 
the same galaxy (UGC 929). In the two sets of plots (a) and b) ), the left panel displays 
all 1080 points and the right panel shows only the (< 50) points remaining after iterative 
filtering. Set a) shows a reasonably well-behaved solution favoring n = 0.6 while set b) shows 
a rather noisy solution favoring n = 0.8. 
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Fig. 8. — Examples of xLner' X^giobai' distributions for a solution with a well-behaved 
Xinner^^ but a flat x^obai' distribution (UGC 784 B-band), plot a), and for a very well-behaved 
solution in both distributions (UGC 929 B-band), plot b). 
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Fig. 9. — Comparison of different bulge fits for tlie same profile (UGC 784 B-band). The 
plot on the right has a bulge fit (dashed-dotted blue line) which is likely unphysical. Its x^u 
however, is lower than that of the decomposition on the left plot, whose bulge fit looks more 
realistic. Without adopting the xL statistic, the plot on the right is favored. Using the Xin 
in addition to the x^i &s a discriminator, the plot on the left is favored. (See left plot of 
Fig. 8 for the corresponding x^ vs. n distributions.) Symbols, colors and line- types are as 
defined in Fig. 11. 
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Fig. 10. — Histograms of sky and seeing FWHM offsets preferred in our analysis for all 
profiles surviving the final cut, separated into the four different bands. Note that the H- 
band sky error is more than an order of magnitude smaller than in the optical (as in the 
actual measurements). 
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Fig. 11. — Decomposition results for a Type I galaxy (UGC 9908). In the upper panels 
of each plot, the data points and measured sky error envelopes are shown with solid black 
circles and dots respectively. The blue dashed-dotted and green dashed lines show the bulge 
and disk fits respectively, and the solid red line is the total (bulge+disk) fit. The fits are 
all seeing-convolved using the best selected seeing values. The bottom panels show the fit 
residuals where A/i(r) = data{r) — model{r). 
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Fig. 12. — Decomposition results for a Type-II/Transition galaxy (UGC 12527). 
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Fig. 13. — Decomposition results for a galaxy with a truncated disk (UGC 927). Note that 
sky errors could not account for the truncation. 
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Fig. 14. — Decomposition results for a galaxy with a "bulgeless" disk (UGC 10757). 
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Fig. 15. — Histograms of Sersic n parameter for "final" solutions (left), and the reduced set 
of solutions after further visual examination (right). See text for details. 
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Fig. 16. — Histograms of Sersic n parameter fitting n as a free parameter in the decompo- 
sitions. Results using three different values for the initial estimate of n are shown: n = 0.2 
(left), n = 1.0 (middle), n = 4.0 (right). Note the different y-axis scales in each of the plots. 
The selection criteria for the fits is as described in the text. 
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Fig. 17. — comparison of floated n versus fixed n solutions. The point types and colors are 
as follows: B-band (triangles), V-band (squares), R-band (pentagons), H-band (asterisks), 
Type-I (blue), Type-II (red), and Transition (green). Note the different axis scales for Xin- 



- 64 - 



2 - 







2 - 



1 1 1 1 1 1 III 1 


1 1 1 1 - 

• 

H Ji 
• — 

• 

» • 

■ ■ - 

> • H - 
^ ■ ■ • 

» • 

> • ▲ ▲ ■ 

• • 

I • ▲ • • - 
» • ■ ■ 

« ▲ A - 
▲ 

▲ ▲ 


1 1 1 1 1 1 1 1 1 1 1 1 1 1 


B - 

A 

» • 

• • • 
» • 

• 

» • ▲ 

t -• • • -• 

» • • • 
> • • 

I • • - 

• « • • 

• A A - 
A ■ 

T 1 t , " 







6 



Morphology 



Fig. 18. — Sersic n versus morphological type index. Blue circles, red triangles, and green 
squares indicate Type-I, Type-II, and Transition galaxies respectively. 
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Fig. 19. — Bulge and disk parameters versus ellipticity (1 — b/a). The point types are as 
follows: Type-I (blue triangles), Type-II (red squares), Transition (green circles). 
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Fig. 20. — versus for our current Type-I data (solid symbols) and the decompositions 
of Graham (2001) of de Jong & van der Kruit (1994) 's data (open symbols). Blue circles 
are B-band, green pentagons (our data only) are V-band, red triangles are R-band, and 
magenta squares are H-band (us) and K-band (Graham (2001)). The dashed lines have a 
slope {velh) = 0.22 for late-type spirals. Note that the large dispersions in the and 
(Table 3) counteract to yield significant Ve/h correlations. The left plot is in apparent units 
(arcsec) and the right plot shows the physical scale in kpc. The discrete nature of our data 
in the right plot is due to the limited precision of the measurement (one decimal). 
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Fig. 21. — Distribution of r^/h with Hubble types for our Type-I galaxies and those of 
Graham (2001). Symbols and colors are as in Fig. 20. The dashed line describes the fit 
(re//i) = 0.20 — 0.013(T — 5) with la = 0.09 errors (dotted lines) based on our data only. 
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Fig. 22. — Distribution of the Sersic n parameter from cosmological simulations by Scanna- 
pieco & Tissera (in prep.). 
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Fig. 23. — Difference between the exact numerical value for 6„ and several commonly adopted 
approximations. The short (red) and long (green) dashed lines are the two most commonly 
used approximations found in the literature. The solid blue line shows Ciotti & Bertin's 
asympotic expansion and the dotted purple line depicts our adopted extension at n < 0.36. 
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Tablc 1. Tabic of mean values and mean rms deviations for repeat observations. 





n 




(mag/D ) 




re 

(kpc) 




(mag/D ) 




h 

(kpc) 




N 


I B 


1.00 


0.13 


21.88 


0.19 


1.13 


0.08 


20.90 


0.03 


5.35 


0.18 


7 


V 


0.98 


0.07 


20.16 


0.27 


0.81 


0.11 


19.96 


0.03 


4.06 


0.08 


7 


R 


0.88 


0.16 


20.02 


0.17 


0.77 


0.11 


19.59 


0.04 


3.76 


0.07 


13 


H 


0.95 


0.14 


17.52 


0.26 


0.85 


0.14 


17.33 


0.10 


3.49 


0.17 


5 


Total 


0.94 


0.13 


20.07 


0.21 


0.87 


0.11 


19.61 


0.05 


4.13 


0.11 


32 


II B 


0.72 


0.17 


22.84 


0.28 


1.05 


0.36 


21.45 


0.06 


3.60 


0.13 


2 


R 


0.61 


0.19 


20.00 


0.34 


0.41 


0.05 


19.80 


0.02 


3.01 


0.02 


4 


H 


0.68 


0.25 


17.61 


0.07 


0.80 


0.00 


17.45 


0.09 


5.28 


0.39 


2 


Total 


0.66 


0.20 


20.12 


0.26 


0.67 


0.12 


19.62 


0.05 


3.73 


0.14 


8 


Tr V 


0.65 


0.14 


21.41 


0.27 


0.78 


0.11 


20.69 


0.05 


5.00 


0.14 


2 


R 


0.65 


0.11 


20.19 


0.33 


0.71 


0.37 


19.60 


0.03 


4.51 


0.06 


5 


H 


1.18 


0.40 


18.14 


0.41 


0.60 


0.09 


17.27 


0.04 


2.58 


0.17 


3 


Total 


0.81 


0.20 


19.82 


0.34 


0.69 


0.07 


19.12 


0.04 


4.03 


0.11 


10 
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Table 2. Comparison of B/D decomposition parameters. 



Author 


UGC 


band 


1^0 


h 




re 


n 








(mag/D ) 


n 


(mag/n ) 


n 




Graham 


463 


B 


20.77 


13.0 


20.51 


1.6 


0.4 


de Jong 


463 


B 


20.76 


13.5 


20.60 


1.3 


1.0 


us 


463 


B 


20.59 


14.2 


20.15 


1.1 


0.4 


de Jong 


463 


TT 

H 


16.80 


12.0 


16.73 


1.9 


1.0 


us 


1(33 


H 


17.11 


13.7 


16.56 


1.6 


0.5 


Graham 


oUoU 


R 
D 


00 on 
zz.zu 


1 Q 
iy.o 


00 AO 
ZZ.'iZ 


1 A 


n Q 
u.o 


de Jong 


3080 


B 


21.99 


17.2 


19.88 


0.2 


1.0 


us 


3080 


B 


21.99 


24.3 


22.74 


2.0 


0.5 


de Jong 


3080 


H 


18.21 


15.1 


18.79 


1.7 


1.0 


us 


3080 


H 


18.28 


18.2 


19.21 


2.7 


0.9 


Graham 


3140 


B 


21.06 


13.7 


22.16 


4.8 


2.9 


de Jong 


3140 


B 


20.90 


13.1 


20.46 


2.0 


1.0 


us 


3140 


B 


20.52 


12.9 


20.75 


2.8 


1.1 


us 


31 10 


B 


20.66 


13.6 


21.28 


3.6 


1.1 


de Jong 


3140 


H 


16.96 


11.3 


15.99 


2.1 


1.0 


us 


3110 


H 


17.22 


12.0 


17.28 


1.1 


1.9 



Note. — Parameter comparison between de Jong (1996a), Graham 
(2001), and us. Graham's decompositions (B-band only as Graham 
did not re-analyse de Jong's H-band data) use de Jong's ID brightness 
profiles. Our decompositions are based upon data discussed in this 
paper and presented in Paper II. The surface brightnesses have not 
been corrected for Galactic or internal extinction in order to compare 
with de Jong (1996a) and Graham (2001). Repeat observations are 
hsted for UGC 3140 B-band. The galaxies UGC 463 and 3080 are 
Type-II systems; UGC 3140 is a Type-I galaxy. See Graham (2001; 
Fig. 6) for examples of his decompositions (esp. for UGC 463). 
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Table 3. Mean and standard deviation, for bulge and disk parameters for different 

types and bandpasses. N denotes the number of data points. Multiple observations count 

as independent entries. 



Par am 


Type 




B 


V 


R 


H 




I 


N 


56 


39 


59 


34 




II 


N 


23 


23 


31 


22 




Tr 


N 


4 


12 


19 


16 


n 


I 


u((t) 


0.88 (0.38) 


1.02 (0.42) 


0.95 (0.35) 


1.11 (0.57) 




II 


fi{a) 


0.63 (0.33) 


0.65 (0.25) 


0.67 (0.23) 


0.68 (0.32) 




Tr 




0.50 (0.37) 


0.89 (0.41) 


0.79 (0.40) 


1.04 (0.50) 




I 




1.00 (0.77) 


1.07 (0.81) 


0.96 (0.62) 


0.87 (0.52) 


(kpc) 


II 




0.72 (0.41) 


0.65 (0.32) 


0.66 (0.32) 


0.65 (0.36) 




Tr 




0.40 (0.14) 


0.68 (0.27) 


0.69 (0.31) 


0.70 (0.33) 


h 


I 




5.48 (2.63) 


5.03 (2.40) 


4.54 (2.08) 


3.81 (1.58) 


(kpc) 


II 




4.84 (1.85) 


4.61 (1.81) 


4.23 (1.65) 


4.17 (1.96) 




Tr 




4.05 (1.07) 


4.51 (1.58) 


4.34 (1.73) 


3.26 (1.25) 


r /h 


I 




IQ fO OQ^I 


21 fn ini 


22 fO 09") 


23 fO 09") 




II 


u((t) 


0.16 (0.10) 


0.14 (0.05) 


0.16 (0.06) 


0.16 (0.08) 




Tr 


ll{ CTx 
) 


10 Co 021 


16 fO 08") 


1 6 fO 04") 


21 fO 07") 




I 




21.76 (0.95) 


20.69 (0.99) 


20.13 (0.94) 


17.51 (1.03) 


(mag/D") 


II 




22.11 (1.04) 


21.12 (0.80) 


20.42 (0.95) 


17.62 (1.16) 




Tr 


H{a) 


22.21 (0.91) 


21.57 (0.82) 


20.61 (0.92) 


18.11 (1.62) 




I 




21.15 (0.65) 


20.26 (0.63) 


19.75 (0.56) 


17.28 (0.57) 


(mag/D") 


II 




21.42 (0.50) 


20.54 (0.45) 


19.98 (0.52) 


17.48 (0.46) 




Tr 


fi{a) 


21.27 (0.69) 


20.45 (0.62) 


19.83 (0.64) 


17.42 (0.91) 




I 


H{a) 




1.07 (0.09) 


1.15 (0.17) 


1.34 (0.21) 






N 




53 


77 


49 




II 


ii{a) 




1.06 (0.07) 


1.15 (0.09) 


1.36 (0.17) 






N 




19 


23 


13 




Tr 


H{a) 




1.01 (0.08) 


1.09 (0.10) 


1.34 (0.09) 






N 




4 


3 


3 
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Table 4. Decomposition Results for Type I Profiles 



Profile 


e 


Sky 


FWHM 


n 


Me 






Mo 


h 


h 


B/D 


UGC obs band 


(l-6/a) 


(mag/n ) 


(") 




(mag/D ) 


(") 


(kpc) 


(mag/D ) 


(") 


(kpc) 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 



10 


044 


B 


0, 


.295 


22.33- 


1, 


.8- 


1, 


n+0.5 

■^— 


91 Q9+0.59 
^^•^^-0.00 


2, 


9+0.6 
' — 0.0 


1, 


7+0-5 
■ ' —0.0 


22.011°°^ 


15, 


9+1.1 
— 0.0 


11.81°'° 


04+0-01 


10 


043 


V 


0, 


.295 


21.55- 


2 


.0 


1, 


4+0-5 

—0.0 




2, 


0+1.6 

■°-o 


2 


9+1.3 

■ — 


21.09l°:J8 


14, 


0+1.6 

—0.0 


11-liao 


0.07l°-°3 


10 


042 


R 


0, 


,295 


21.09- 


1 


.8 


1, 


p:+0'5 
"-^—0.0 


20.82l°;80 


3, 


0+2.1 
'"—0.0 


2 


4+1-7 
■ — 0.0 


20.50t°-§g 
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7+l'6 

' —0.0 


lO'Slao 


0-081°°^ 


10 


754 


H 


0, 


.295 


13.96- 


1, 


.2 


1, 


q+0.3 
"^-0.1 


17.48t°-j8 


1, 


7+0.5 
' ' —0.1 


1 


0+0.4 

■'-^ — 0.1 


17.58i;;;?« 


9, 


9+0.7 
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7 1+0.5 
' ■ — 0.0 


0.08l°'°J 


16 


036 


B 


0, 


.157 


22.25+ 


1, 


.4- 


1, 


9+0.1 
'^ — 0.1 


20, 


qn+o-09 

'^"-0.15 


2, 


q+0.0 
'''—0.2 


1 


0+0.0 

■"-0.1 


21 36+"'"" 

Z1.0U_Q 03 
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4+0.0 

'^ — 0.5 


"■^—0.2 


0.08l'>:°° 


16 
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V 


0, 


,157 


21.44+ 


2 


.1 


1, 


9+0-1 

— 1 


21, 


ns+0-15 


6, 


0+0.6 
'°-0.6 


2 


4+0-2 
■ — 0.2 




20, 


q+0.1 
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7 -1+0.0 
' ■ — 2 


f^ 99 + 0.02 

^J'-^-^ — O 02 


16 


034 


R 


0, 


.157 


20.53 


1, 


.4- 


1, 


4+0-0 

— 0.2 


19, 


.c:+0.00 
'^^-0.31 


3, 


9 + 0.0 
' — 0.3 


1 


1 +0.0 
■-"--0.1 


19.88«'«? 


16, 


c+O.O 
"-^-0.8 


5.8l°'0 


12+°'°l 


52 


026 


B 


0. 


,175 


21.93 


2, 


.1 


0, 


7+0.1 
— 0.0 


22, 


it;+0.14 
'l^'-O.OO 


2, 
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'"—0.0 





7+0.0 
■ ' —0.0 


21.10«'«1 


13, 


q + 0.3 

'-'-0.1 


4.7lg:S 


01+"'°" 
"•"1 — 0.00 


52 
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V 


0. 


,175 


21.50- 


1, 


.2 


0, 


O+O-l 

■°-o.o 


21, 


9f-+0.12 
'^"-0.00 


3, 
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'"— 


1 
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20.16l°-°« 


12, 
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'" — 


4 2+°-l 
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52 


025 


R 


0, 


,175 


20.34- 


2 


.3+ 


1, 


1+0.1 

— 1 
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'="^-0.12 


3, 
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1 
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— 
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^■^^ — 0.00 


52 
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R 


0, 


,175 


21.07 
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.8 


1, 
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2, 
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1, 


0+0.1 


19.57tHg 


11, 


4+0-2 
—0.0 


4 0+°-i 
^•^J— 0.0 


04+° °l 
"•u^— 0.00 
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031 


B 


0, 


,197 


22.05+ 


2 


.9+ 


0, 
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■^—0.1 


20.05l°:?| 


1, 
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14 1+!!-;! 


3 2+!!-! 
"•^-0.0 


03+::-?! 
"•""-0.00 


3062 


089 


R 


0.487 


21.04 


1.6- 


1 


•O+nn 
— 0.{) 


21.72 + nni 
' —0.00 


5, 


-1+0.7 
' — 0.0 


1. 


9+0.2 
■^-0.0 


19 •59+0-03 

— 0.00 


13.4+!!-? 

— 0.0 


3.1+!!-!! 

— 0.0 


0.03+°-?! 

— 0.00 


3062 


534 


H 


0.487 


13.61- 


0.9- 





^ — 0.0 


17.52+°-3,2 

' ^ — 0.06 


1, 


4+0-1 
• — 0.0 


0. 


Q + 0.0 

■•^-0.0 


16.64+°-°! 

^ " — 0.00 


9.8+°-! 

" — 0.0 


9 9+0.0 
-0.0 


0.01+2-?? 

— 0.00 


3070 


094 


B 


0.397 


22.45+ 


1.9- 





7+0.1 

■ —0.0 


22.79+°7,n 

— 0.00 


5, 


2+0.2 

' —0.2 


0, 


0+0.0 
■8-0.0 


21.41+!!-!!! 

— 0.02 


12.6+°-!! 

— 0.2 


2.0+°-? 

— 0.0 


0.08+?-?° 

— 0.00 


3070 


093 


V 


0.397 


21.35+ 


2.1j- 

-r 


0, 


^—0.1 


22.20'^n'?n 
■^■^■■^'-'—0.10 


5, 


c+0.0 
■"—0.5 


0. 


n+0.0 
■9-0.1 


20.95+!!-°° 


12.5+°-!! 

-^■^•'-'—0.6 


2.0+!!? 

^—0.1 


0.11+°-°? 

^ — 0.01 


3070 


092 


R 


0.397 


20.81+ 


1.8- 





o+O.l 
•°— 0.0 


21.58+°?,2 

^"—0.00 


5, 


9+0.4 
• — 0.2 


0. 


Q+0.1 
■8-0.0 


20.57+°-°t 


12.2+°-2 
^•'■^— 0.6 


2.0+°-? 

^—0.1 


0.12+°-0J 
'-'•-^'^—0.01 


3140 


079 


B 


0.076 


20.95- 


1.7- 


1 


(,+0.2 
• — 0.0 


20.66'^nnn 
^^^ — 0.00 


2, 


7+0.4 
' —0.0 


0, 


0+0.1 
■«-o.o 


20.45'^nnn 

^ ^ — 0.00 


11.5+!-? 


3.4+?-^ 

^ — 0.0 


0.09+°-°! 

^^^ — 0.00 


3140 


066 


B 


0.076 


22.07- 


1.7+ 


1 


1+0.1 
•-*- — 0.1 


20.71+°-;,^ 

^ ' —0.09 


2, 


0+0.4 

■"—0.0 


0, 


0+0.1 
■«-o.o 


20.48"^nnn 
^ " — 0.00 


"^ — 0.0 


3.8+!!-! 

" — 0.0 


0.08+°-°! 

^^" — 0.00 


3140 


061 


B 


0.076 


21.90 


2.5+ 


1 


•^-0.1 


21 24+°-?5 
^-^•^^-0.15 


3, 


6+°-? 
"-0.3 


1, 


-1+0.1 
■1-0.1 


20 62+!:-!l;! 
^"•"^-0.03 


13 6+!!-? 

-'-"•"-0.3 


4 0+!!-! 
^•"-0.1 


09+?-?! 

"•""-0.01 


3140 


096 


V 


0.076 


20.64 


2.6 


1 


c+O.O 
• — 0.3 


20.34+°-°° 

— 0.36 


4, 


0+0.0 
■ — 0.7 


1, 


q+0.0 

■•^-0.2 


20.04+!!-!!!! 

— 0.09 


13.9+''-? 

^"^ — 0.6 


4.1+''-! 

— 0.2 


0.17+°-?° 

^ ' —0.02 


3140 


060 


V 


0.076 


20.88 


2.4 


1, 


4+0.1 
•^—0.1 


20.18+°-l;i 


4, 


-o+S-? 

^—0.2 


1, 


9+0.1 
■^-0.1 


— 0.03 


13.9+!!-* 

" — 0.4 


4.1+!!-! 

— 0.1 


0.16+°-°! 

^ ^ — 0.00 


3140 


078 


R 


0.076 


20.44 


1.8- 


1, 


O+0.2 
— 0.0 


19.37+°-24 

— 0.00 


3, 


g+0.6 
—0.0 


1, 


1+0.2 
■1-0.0 


19.40+!!-" 

— 0.00 


12.8+!!-? 

— 0.0 


0+0.2 
— 0.0 


0.17+?-?? 

— 0.00 


3140 


065 


R 


0.076 


20.66- 


1.8- 


1 


O+0.2 
■'^ — 0.0 


19.36+°-2* 


3, 


C.+0.5 
—0.0 


1, 


-1+0.2 
■1-0.0 


19.35+°-°^ 

^^^ — 0.00 


12.6+°-6 

^ — 0.0 


3.7+°-? 

^ —0.0 


0.16+°-°? 

^ ^ — 0.00 


3140 


059 


R 


0.076 


20.30 


1.8- 


1 


O+O.l 

— 0.0 


19-43+°-l,2 

— 0.04 


3, 


,6+°-? 

—0.1 


1, 


1+0.1 

■-'--o.o 


19.42+°-!!? 

— 0.05 


13.0+°-! 

— 0.5 


3.9+°-? 

— 0.1 


0.16+::-?! 

— 0.00 


3140 


448 


H 


0.076 


14.31+ 


2.2+ 


1, 


.9+0-1 

^—0.1 


17.19+°-?* 

— 0. 10 


4, 


1+0.9 
—0.1 


1, 


9+0.3 
■^-0.0 


17.13+!!-!* 

— 0.03 


12.0+°-? 

^■^•"—0.3 


3.6+°-! 

"^—0.1 


0.28+°-°? 

" ^ — 0.00 


3245 


042 


B 


0.287 


21.73+ 


1.6- 





O + O.l 

■"—0.0 


21.37+°-2g 

' — 0.08 


1, 


0+0.0 

■"—0.2 


0. 


(5+0.0 
■"-0.1 


21.65+°-°° 

— 0.02 


18.0+°-° 

— 0.5 


5.8+0-0 

— 0.2 


0.02+0-00 

" " — 0.00 


3245 


041 


R 


0.287 


20.25+ 


2.5+ 


1 


.o+°-S 

— 0.2 


20.03+°-?9 

— 0.14 


1 


^+0.1 
—0.1 


0, 


c+O.O 

■°-o.o 


20.24+°-°° 

— 0.02 


15.5_o g 


5.0+0-? 

— 0.3 


0.02+?-?? 

— 0.00 


3245 


030 


R 


0.287 


20.37+ 


1.7- 


1, 


2+0.0 
• — 0.2 


20.52+°-°° 

" " — 0.32 


2, 


1+0.0 
—0.3 


0, 


-+0.0 
■'-0.1 


20.28+°-°° 

" ^ — 0.06 


15.7+°-° 

' — 1.0 


5.1+0-0 

" — 0.3 


0.03+0-00 

^•^^ — 0.00 


3245 


042 


H 


0.287 


13.97+ 


1.8+ 





o+O.O 
•°— 0.1 


17.07+°-°! 


1, 


.4+°-° 

■^—0.1 


4+°-° 
"■^-0.0 


17.70+°-°° 


11.6+°-° 

^—0.4 


3.8+°-° 

" "— 0.1 


0.04+0-00 

^ " — 0.00 


3379 


112 


B 


0.523 


20.75+ 


1.7- 


1 


1+0.2 
•-^-0.0 


20 60+5!-?? 

^"•""-0.08 


1 


g + O.l 

■"—0.1 


0, 


c+O.O 
■^-0.0 


21 06+°-°! 

^-'-•""-0.00 


21 5+°-! 

^-'-•'^-0.1 


6 0+°-? 

""-0.0 


02+?-?° 

"•"^-0.00 


3379 


111 


V 


0.523 


20.47 


1.4 


1, 


2+0.1 
• — 0.1 


19.20+°-" 

— 0.09 


1, 


■8+^-i 

—0.0 


0, 


rr+0.0 

■"-0.0 


19.95+°-°? 


17.6+°-5 

' —0.0 


4.9+°-! 

^ — 0.0 


0.04+°-?° 

— 0.00 


3379 


110 




0.523 


20.05+ 


1.9+ 


1, 


4+0-0 

■ —0.2 


18.96+°-°° 

^ — 0.25 


2, 


,0+°-° 

—0.2 


0, 


g+0.0 


19.46+°-°° 

— 0.03 


18.0+°-° 

— 0.4 


5.0+°-? 

— 0.1 


0.04+°-°° 

— 0.00 


3471 


044 


B 


0.375 


21.77- 


2.2+ 


1, 


t;+0.2 
• — 0.1 


22.64+S'?« 
^ — 0.16 


5, 


7+2.0 
■ —0.6 


2, 


-1+0.7 
■1-0.2 


21.51+!!-!? 

— 0.03 


19.4+!-^ 

— 0.0 


7.2+°-^ 

— 0.0 


0.07+?-?? 

— 0.01 


3471 


033 


B 


0.375 


22.15+ 


2.7+ 


1 


5+0.2 
—0.0 


22.80+ nnn 


6, 


0+2.4 
'—().() 


2, 


q+0.9 

■•^-0.0 


21.52+!!-!? 

— (J. 00 


19.6 + nn 
— (J.O 


7.3+°-* 

— (J.O 


0.07+:;-?* 

—(J. 00 


3471 


804 


H 


0.375 


13.60- 


1.2+ 


1 


1+0.0 
• ^-0.2 


17 7.1+0.00 
-^-'•''^-0.34 


2, 


g+0.0 

■"-0.fi 


1, 


g + 0.0 
■"-0.2 


17 Qi+O.OO 
1 ' •■-'^-0.12 


-'--'-•''-0.6 


4 ■i+o-i 
^••^-0.2 


o.o8t:;:;;° 


4227 


065 


B 


0.115 


22.10- 


1.3- 





5 + 0.1 

•°-o.o 


20 93+0-1" 


3, 


9 + 0.0 
■^-0.1 


0, 


+ 0.0 

■*-o.o 


9-1 qq + ().03 


14 1+0-5 
^^•1-0.0 


Q g + O.l 

"•°-o.o 


12+"-°° 

'^•1^-0. 00 


4227 


064 


R 


0.115 


20.89+ 


1.9- 


0, 


O+O.l 

•^-0.0 


19 92+°-°* 


3, 


g+0.0 

■°-0.2 


0, 


g+0.0 
■"-0.1 


20 26+° °^ 


10.C)_g 7 


Q q+0.1 

0.»_g 2 


Q 1Q+0.00 

>J^1"-0.01 


4227 


021 


H 


0.115 


14.20- 


2.8+ 


0, 


Q + 0.1 

•^-0.0 


17 77+0.17 
1'-' '-0.00 


3, 


g+0.4 

■°-o.o 


0, 


g + O.l 

■^-0.0 


i8.i5i°;JJ 


l3.9l!;J 


Q K+o.a 
"•"-0.0 


17+001 

U-l'-O.OO 


4227 


123 


H 


0.115 


13.93- 


3.2+ 





q+0.1 
•^-0.0 


17.83l°:J« 


3, 


(5+0.6 
■''-0.0 


0. 


q+0.2 
■^-0.0 


i8.i3l°:Jg 


13.8lJ:g 


Q c+O.S 
"■"-0.0 


-,g+0.02 

<-'.lD_g gg 


4778 


071 


B 


0.086 


19.39+ 


0.3 





(5+0.0 
•O-O.O 


20.781°-°^ 


2, 


(5+0.0 

■''-o.o 


1, 


q+0.0 
■^-0.0 


91 99+0.00 

Z1.ZZ_Q gg 


ii.ot°-° 


8.0t°-°2 


lS+0 02 
"•I'-'-O.OO 


4778 


070 


R 


0.086 


19.26+ 


2.3+ 


1, 


t;+0.1 
•°-0.2 


20-19l°-i«8 


5, 


o + l-l 
■■^-1.6 


3, 


q+0.8 
■"-1.2 


20 si+°-'!° 


11.2t°-! 


8.il^:J 


57+0-26 
^•0' -0.24 


4978 


090 


B 


0.225 


22.08+ 


1.5- 


1 


Q + 0.0 

•'-'-0.0 


23.77tO-°,f, 


8, 


t; + 0.0 
■"-0.6 


2, 


9+0.0 
■^-0.2 


9q p.q+0.00 


23.5+°-; 


6.0l?-° 


20+0-02 


4978 


089 


R 


0.225 


20.70- 


2.3+ 





O+O.l 

•°-o.i 


22 13+0-13 


2, 


0+0.6 
■=-0.2 


0, 


7+0.2 
■'-0.0 


2l.29t»;;« 


10.2t°i 


2 g+0.2 
^•"-0.1 


o.o6t:;-°,i 


5674 


092 


B 


0.296 


21.98 


1.8- 


1, 


Q+0.2 
••^-0.0 


,1 Q2+0.30 
^^•^^-0.00 


3, 


4+0.9 
■^-0.0 


2, 


1+0.6 
■1-0.0 


21 4t;+0-08 

Z1.^0_g QQ 


i4.5l°;J 


Q+0.4 
8-9-0.0 


o.08l°o;°J 



-75- 
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Profile 


£ 


Sky 


FWHM 


n 


Me 


re 


re 


MO 


h 




B/D 


UGC obs band 


(1-6/a) 


(mag/D ) 


(") 




(mag/D ) 


(") 


(kpc) 


(mag/D ) 


(") 


(kpc) 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 



5674 


054 


B 





296 


22.06+ 


2.0- 


1 


,+0.3 
"^-0.0 


22 


nQ+o.3s 

^■^-0.00 


3 


4+0.9 
4-0.0 


2 


1+0.6 
1-0.0 


91 ,0+0.06 
Zi.OO_o Q2 


13 


q+0.3 

^-0.2 


8 


5+0.2 
0-0.1 


07+°-°2 
'3-l"-0.00 


5674 


091 


R 





296 


20.57- 


1.8- 


1 


4+0.4 


20 


,9+0.60 
•^^-0.00 


3 


0+2.3 

«-0.0 


2 


, + 1.4 

•'-o.o 


20 2f\+^-'^'^ 

ZU.ZD_(, QQ 


12 


g+3.2 
'3-0.0 


7 


7+2.0 
'-0.0 


19+0-13 

"J-i^-o.oo 


5808 


078 


B 





085 


22 


14+ 


1.8- 


1 


Q + O.fi 


22 


01+0.84 
^^-0.00 


1 


7+1.0 
' -0.0 





q+0.5 
"-0.0 


20 93+0-03 
^"■■'■'-0.00 


8 


5+0.1 
^-0.0 


4 


9+0.0 
^-0.0 


01+° °^ 
"■'31-0.00 


6128 


085 


B 





164 


21 


80 


0.4+ 





7+0.1 
'-0.0 


20 


c:9+0.13 
■-'^-0.00 


1 


6+0.1 
"-0.0 





1 +0.0 
1-0.0 


20 21+"-"^ 


12 


5+0.2 
^-0.0 


1 


1+0.0 

1-0.0 


02+° °° 
"■"3^-0.00 


6128 


084 


R 





164 


20 


27+ 


1.2- 


1 


,+0.2 
■^-0.0 


19 


69+0-28 
"^-0.04 


2 


C.+0.4 
^-0.1 





9+0.0 
^-0.0 


18 96+''-"2 
io.yD_(, g2 


11 


4+0-1 

4-0.0 


1 


Q + 0.0 


0-05l°:SJ 


6128 


235 


H 





164 


13 


72+ 


1.4- 


1 


O+0.2 
^-0.1 


17 


84-0. 19 


3 


,, + 1.2 
''-0.5 





-1 + 0.1 
'^-0.0 


1 e. -71 +0.11 

16-71-0.06 


10 


7+0.4 
'-0.3 





q+0.0 
"-0.0 


0.08l'>:- 


6413 


080 


B 





125 


21 


52+ 


2.3- 


1 


(,+0.2 


22 


S/1 + 0-22 


3 


g+0.3 
"-0.2 


1 


c+O.l 
■^-0.1 


21.67±ng 


11 


0+0.1 


4 


7+0.0 
'-0.3 


n 07+°-°i 

"■'3'-0.00 


6413 


079 


R 





125 


20 


00+ 


2.1- 


1 


1+0.1 
^-0.1 


21 


it:+0.13 
^■^-0.09 


3 


0+0.4 
^-0.5 


1 


g+0.2 
"-0.2 


20.65in? 


10 


5+0.1 
0-1.1 


4 


5+0.0 
0-0.5 


16+°°3 
UlD-0.02 


7357 


089 


B 





251 


21 


70- 


2.4+ 





g+0.3 
°-0.0 


22 


44+0-34 

^^-0.00 


2 


9+1-7 

'^-O.O 


1 


9+0.7 
''-O.O 


21 74+0-30 

^1- '4-0.00 


13 


9+2.5 
^-0.0 


5 


g+1.1 
°-0.0 


04+0 03 

"•'34-0.00 


7357 


088 


R 





251 


20 


23 


1.7- 





c+0.1 
^-0.0 


20 


qc+O-OS 
'^'-'-0.00 


2 


q+0.2 
^-0.1 


1 


9+0.1 
^-0.0 


90 56+°'l° 

ZU.OD_g gg 


10 


7+0.5 
' -0.2 


4 


5+0.2 
O-O.l 


Q Q~+0.01 

"■"' -0.00 


7357 


077 


R 





251 


20 


69- 


2.0+ 





fj+0.1 


20 


04+0.13 

°4-0.00 


2 


t;+0.4 

■'-O.Q 


1 


(,+0.2 


20.58j:g:oo 


10 


7+0.9 
'-0.0 


4 


5+0.4 
O-O.O 


07+° °l 
"■'3'-0.00 


7632 


036 


B 





431 


22 


11+ 


0.4+ 





Q + 0.0 

^-0.1 


21 


t;fi+0.00 
■^°-0.12 


2 


,+0.0 
•^-0.2 


1 


1 +0.0 
-■--0.1 


21.861-;" 


18 


0+0.0 
^5-0.6 


8 


q+0.0 


04+° °° 
"■'34-0.00 


7632 


035 


V 





431 


21 


21+ 


1.9+ 


2 


(,+0.1 


21 


94+0-18 
^4-0.19 


4 


4+0.6 
4-0.6 


2 


1 +0.3 
1-0.3 


2i.i6l!!:r4 


18 


g+0.1 
'3-0.8 


8 


0+0.0 
8-0.4 


07+° °l 
"■'3' -0.00 


7632 


034 


R 





431 


20 


73+ 


1.7+ 


1 


g+0.2 
°-0.1 


20 


79+0.32 
'^-0.15 


3 


9+0.7 
^-0.3 


1 


c+().3 
^ -0.1 


2o.55i°:;;| 


16 


g+0.0 
"-0.7 


8 


Q + 0.0 

^-0.3 


07+°-°i 

"■"' -0.00 


7632 


223 


H 





431 


14 


48- 


2.7+ 


1 


a+0.3 
**-0.1 


18 


o,+0.53 

^■^-o.ie 


4 


1 +2.0 
1-0.2 


1 


g + 0.9 

•^-0.1 


18-12li',;oo 


13 


9+1.1 
^-0.0 


6 


9+0.5 
^-0.0 


l2+»-°5 
"■1^-0.01 


9467 067 


B 





431 


22 


53+ 


1.6 





t; + 0.1 
^-0.0 


22 


ro+0.14 
3*^-0.03 


2 


1 +0.2 
-^-0.1 





c+O.O 
^-0.0 


2o.9oi;j:»i 


11 


9+0.0 
^-0.1 


2 


7+0.0 
'-0.0 


01+° °° 
"■"1-0.00 


9467 


068 


B 





431 


22 


53+ 


1.6 





c+O.l 
^-0.0 


22 


CQ + 0.14 


2 


1 +0.2 
-^-0.1 





c+O.O 

''-o.o 


2o.9oi;j:»i 


11 


9+0.0 
^-0.1 


2 


7+0.0 
' -0.0 


01+° °° 

"■"1-0.00 


9467 


064 


B 





431 


22 


46+ 


2.0+ 





o+O.l 
°-0.2 


21 


70+0.13 
'°-0.14 


1 


c:+0.2 
'^-O.O 
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